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ABSTRACT 


The  past  twenty  years  have  seen  an  explosion  in  the  realm  of  infrared  deteetion 
teehnology  fueled  by  improvements  in  III-V  semieonduetor  teehnology  and  by  new 
semieonduetor  growth  methods.  One  of  the  fastest  growing  areas  of  this  researeh 
involves  the  use  of  bandgap  engineering  in  order  to  ereate  artifieial  quantum  wells  for  use 
in  Quantum  Well  Infrared  Photodeteetors  (QWIPs).  QWIPs  have  an  advantage  over 
other  infrared  deteetors  sueh  as  Mereury  Cadmium  Telluride  (MCT)  beeause  they  have 
larger  bandgaps  and  are  therefore  stronger  and  eheaper  to  manufaeture.  This  thesis 
introduees  one  method  of  “multi-eolor”  detection  through  the  use  of  an  asymmetric 
quantum  well  structure  in  which  all  energy  transitions  are  possible.  The  QWIP  structure 
in  this  thesis  was  designed  to  detect  a  laser  wavelength  of  1.06  |im  and  a  wavelength  in 
the  8-10  |im  atmospheric  window. 

The  relevance  of  a  detector  that  is  tuned  to  these  wavelengths  is  that  it  can  be  used 
on  military  aircraft  as  a  laser  spot  tracker  and  an  infrared  imager  providing  much  greater 
accuracy  and  dependability  than  older  systems. 
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I.  INTRODUCTION 


A.  QUANTUM  WELL  INFRARED  PHOTODETECTOR  (QWIP) 

The  past  twenty  years  have  seen  an  explosion  in  the  realm  of  infrared  detection 
technology  fueled  by  improvements  in  III-V  compound  semiconductor  technology  and  by 
new  semiconductor  growth  methods  such  as  molecular  beam  epitaxy  (MBE)  and  metal 
organic  chemical  vapor  deposition  (MOCVD).  These  technologies  have  made  it  possible 
to  create  artificially  structured  superlattices  for  use  in  novel  device  applications.  The 
QWIPs  operate  by  absorbing  incoming  photons  and  exciting  carriers  from  the  ground 
state  to  an  excited  state,  which  is  known  as  intersubband  transitions  and  is  shown 
schematically  in  Figure  1.1.  These  excited  carriers,  or  photoelectrons,  can  then  be  used 
to  create  a  photocurrent  through  the  application  of  an  external  bias. 


Figure  1 . 1  Diagram  showing  quantum  well  structure,  energy  levels energy 

gap  and  photon  absorption  {fico). 

QWIPs  have  an  advantage  over  other  infrared  detecting  materials  such  as 

HgCdTe  (MCT),  which  is  shown  in  Figure  1.2,  because  they  have  much  larger  bandgaps 

and  therefore  are  stronger  and  easier  to  manufacture.  QWIPs  also  have  an  advantage 
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over  MCT  devices  due  to  the  fact  that  MCT  devices  utilize  interband  transitions,  which 
limit  the  ability  of  the  detector  to  be  “tuned”  to  a  specific  wavelength.  These  advantages 
have  made  QWIPs  the  focus  of  recent  infrared  detection  research. 


Figure  1 .2  Diagram  of  HgCdTe  bandstructure  showing  the  energy  gap  and  the  photon 

absorption  across  the  gap. 

The  idea  of  using  semiconductors  to  create  quantum  wells  was  first  proposed  by 
Esaki  and  Tsu  (1969)  and  since  then  many  studies  have  been  done  on  artificial  quantum 
wells,  which  has  led  to  a  detailed  understanding  of  energy-level  spacing,  intersubband 
and  interband  selection  rules,  line  shapes  and  intersubband  absorption.  The  earliest 
working  QWIP  was  designed  using  GaAs/AlGaAs  and  was  based  on  bound-to-bound 
transitions  (Levine  et  al.  1987).  The  responsivity  of  this  first  device  was  extremely  low 
because  the  photoelectrons  had  to  tunnel  through  large  barriers  and  therefore  the  device 
had  to  operate  under  a  CO2  laser.  This  low  responsivity  led  to  the  use  of  bound-to- 
continuum  transitions,  which  was  first  proposed  by  Coon  et  al.  (1984).  The  first  bound- 
to-continuum  QWIP  was  demonstrated  to  have  a  much  higher  responsivity  and  a  larger 
spectral  width  (Hasnain  et  al,  1989),  but  had  a  small  peak  absorbance  due  to  the 
continuum  of  states  above  the  well. 
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Further  advances  in  QWIP  design  have  led  to  structures  that  are  based  on  bound- 
to-quasi-continuum  transitions  (Levine  et  al,  1991),  bound-to-miniband  transitions  (Yu 
et  al,  1991)  and  bound- to-quasi-bound  transitions  (Gunapala  et  al,  1996).  These  more 
recent  configurations  have  shown  good  detector  performance  due  to  the  fact  that  in  all  of 
these  cases  the  ground  state  electrons  do  not  flow  in  response  to  an  external  bias  while 
the  photoelectrons  can  create  a  photocurrent  with  the  use  of  a  relatively  small  bias. 
Multiple  structural  optimizations  have  also  been  made  that  improve  QWIP  performance. 
Increasing  the  barrier  size  has  reduced  the  tunneling  current  by  many  orders  of  magnitude 
(Levine  et  al,  1991),  lowering  the  excited  state  from  the  continuum  into  the  quasi-bound 
region  has  been  shown  to  reduce  the  dark  current  from  thermionic  emission  by  a  factor  of 
~12  at  70K  (Gunapala  et  al,  1996)  and  adding  a  grating  on  top  of  the  device  to  increase 
the  electric  field  polarization  normal  to  the  quantum  wells  in  order  to  increase  the 
absorption  strength  (Karunasiri  et  al,  1995).  All  of  these  QWIPs  have  been  shown  to 
work  well  as  infrared  detectors,  but  they  are  limited  to  the  detection  of  a  narrow  band  of 
wavelengths.  This  limitation  in  QWIP  technology  has  recently  been  overcome  through 
the  design  of  “multicolor”,  or  multiple -wavelength,  QWIPs. 

Multicolor  QWIPs  are  important  because  the  peak  wavelength  emitted  by  an 
object  shifts  as  its  temperature  changes.  This  shift  in  wavelength  can  cause  performance 
limitations  in  single  color  QWIPs  due  to  the  peak  radiation  wavelength  being  far 
removed  from  the  peak  detection  wavelength  of  the  device,  which  will  limit  the 
discrimination  ability  of  the  detector.  One  way  in  which  to  improve  this  limitation  is  to 
design  the  QWIP  to  be  sensitive  to  several  wavelengths  so  that  the  emissivity  changes  of 
an  object  can  be  detected. 

Some  of  the  earliest  methods  of  creating  multicolor  QWIPs  were  to  build  several 

stacks  of  square  quantum  wells,  each  with  a  different  peak  detection  wavelength  (Kock  et 

al,  1992,  Grave  et  al,  1992,  Liu  et  al,  1993).  The  multicolor  response  of  these  stacks 

was  achieved  by  attaching  a  contact  to  each  stack  (Kock  et  al,  1992)  or  by  controlling 

the  external  bias  across  the  device  to  activate  each  stack  sequentially  (Grave  et  al,  1992, 

Liu  et  al,  1993).  Another  method  of  creating  a  multicolor  device  was  to  use  bandgap 

engineering  technology  to  create  an  asymmetric  quantum  well,  where  the  transitions  from 

the  ground  state  to  several  excited  states  are  allowed  and  occur  at  differing  wavelengths 
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(Berger  et  al.,  1992,  Martinet  et  al,  1992,  and  Chiang  et  al,  1996).  The  main  advantage 
of  using  this  system  is  that  it  requires  only  one  set  of  quantum  wells  to  aehieve  multieolor 
detection,  but  the  downside  is  that  in  order  to  get  the  photoelectrons  from  the  lower 
excited  states  out  of  the  well  a  large  bias  has  to  be  applied  to  the  device,  which  increases 
the  tunneling  current.  One  method  of  overcoming  this  large  tunneling  current  is  to  create 
a  superlattice  in  which  an  asymmetric  well  is  surrounded  by  multiple  square  wells  so  that 
the  first  miniband  is  aligned  with  the  lowest  excited  energy  level  in  the  asymmetric  well 
(Ting  et  al,  1997).  The  advantage  of  this  system  is  that  photoelectrons  in  the  lowest 
excited  energy  level  have  a  higher  mobility  due  to  the  presence  of  the  miniband,  and 
therefore  the  responsivity  of  this  structure  is  greatly  improved,  but  the  complexity  of  the 
design  makes  fabrication  more  difficult. 


B,  PURPOSE  OF  THIS  THESIS 

The  purpose  of  this  thesis  is  to  design  an  asymmetric  quantum  well  capable  of 
detecting  a  laser  wavelength  of  1.06  |im  and  a  middle  infrared  wavelength  in  the  8-10  |im 
range  through  the  use  of  both  intersubband  and  interband  transitions.  The  scope  of  this 
research  will  cover  several  aspects  of  semiconductor  technology  including  theoretical 
study  of  band  structure,  intersubband  and  interband  transition  selection  rules, 
intersubband  absorption,  structural  design  of  the  quantum  well,  and  modeling  of  the 
effects  of  an  external  bias  on  the  quantum  well. 

The  solution  to  the  Schrodinger  equation  in  a  QWIP  with  finite  potential  heights 
cannot  be  solved  analytically,  but  must  be  calculated  using  techniques  such  as  the  transfer 
matrix  method  (Vassel  et  al,  1983,  Lui  et  al,  1986,  Wang  et  al,  1989).  The  transfer 
matrix  method  utilizes  the  continuity  and  boundary  conditions  of  the  wavefunctions  to 
obtain  the  energy  eigenvalues  and  wavefunctions  numerically.  These  values  are  then 
used  to  determine  the  absorption  spectrum  of  a  device,  which  is  useful  in  optimizing  the 
device  performance. 

In  order  to  achieve  two-color  detection,  it  is  necessary  to  design  a  quantum  well 
structure  having  two  energy  transitions  corresponding  to  the  two  detection  wavelengths. 

In  a  square  quantum  well,  this  is  difficult  because  intersubband  transitions  can  only  occur 
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between  states  with  opposite  parity  aceording  to  the  selection  rules.  In  addition,  the 
transition  probability  drops  rapidly  as  the  quantum  number  of  the  excited  state  increases. 
Introducing  an  asymmetry  into  the  quantum  structure  and  thereby  altering  the  parity  of 
the  wavefunctions  and  allowing  all  transitions  to  occur  can  overcome  this  restriction 
(Martinet  et  al,  1992).  The  quantum  well  structure  designed  in  this  thesis  makes  use  of 
both  intersubband  and  interband  transitions  to  achieve  multicolor  detection.  This  method 
differs  from  those  studied  in  the  past,  which  utilized  only  intersubband  transitions  to 
obtain  multicolor  detection  (Kock  et  al,  1992,  Grave  et  al,  1992,  Liu  et  al,  1993,  Ting 
et  al.  1997).  The  advantages  of  the  method  used  in  this  thesis  are  that  the  use  of 
interband  transitions  allows  for  much  shorter  wavelength  detection  while  the  longer 
wavelengths  can  be  detected  using  intersubband  transitions.  Furthermore,  the  asymmetry 
of  the  quantum  well  allows  the  excitation  of  electrons  from  the  ground  state  to  the  first 
excited  state  in  the  conduction  band.  This  allows  for  the  extraction  of  the  photoexcited 
electrons  to  occur  at  a  lower  bias,  which  reduces  the  dark  current  dramatically. 


C.  MILITARY  RELEVANCE 

A  multicolor  QWIP  is  relevant  in  the  military  today  because  it  can  be  used  as  both 
an  infrared  imager  and  a  laser  spot  tracker  for  use  in  laser-guided  weapons  delivery.  The 
quantum  well  design  utilized  in  this  thesis  is  capable  of  detecting  wavelengths  as  low  as  1 
|j.m  and  could  easily  be  tuned  for  use  with  any  combat  laser  designator  wavelengths.  In 
addition  to  the  1.06  |im  laser  wavelength,  this  quantum  well  is  designed  to 
simultaneously  detect  infrared  radiation  in  the  8-10  |im  range,  supporting  infrared 
imagery. 

The  applications  of  infrared  and  laser  detection  technology  in  the  military  are 
legion.  Specific  benefits  of  QWIP  detector  devices  are  immediately  evident  in  the 
following  case:  Military  airborne  weapon  systems  often  utilize  a  FLIR  (Forward 
Looking  InfraRed)  for  target  acquisition  and  tracking  and  a  laser  target  designator/range 
finder  that  are  often  co-located  in  a  single  weapons  pod  (e.g.,  LANTIRN)  (Hewish  et  al, 
2002).  Some  aircraft,  such  as  the  F-18,  also  employ  a  laser  spot  tracker,  which  is  located 
in  a  separate,  externally  mounted  pod.  The  addition  of  the  separate  laser  spot  tracker 
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(LSI)  provides  the  pilot  with  a  visual  laser  designation  eue,  whieh  is  espeeially  useful 
when  designation  oeeurs  from  an  offboard  souree  (e.g.,  Speeial  Ops,  “Buddy  lasing”, 
ete.).  The  LST  affords  much  greater  accuracy  and  confidence  in  ordinance  delivery 
while  expanding  employment  options.  However,  boresight  errors  are  inherent  in  a  multi¬ 
pod  system  and  degradations  due  to  misalignment  are  common.  Another  disadvantage  is 
the  loss  of  an  external  weapons  station  for  carriage  of  a  sensor. 

The  quantum  well  detector  design  proposed  in  this  thesis  will  provide  an 
advantage  to  both  airborne  systems,  which  will  lead  to  better  accuracy  and  lethality.  In  a 
system  such  as  LANTIRN,  it  will  provide  laser  spot  tracking  capability  where  one 
currently  does  not  exist.  For  aircraft  such  as  the  F-18,  the  introduction  of  a  laser  spot 
tracker  into  the  weapons  pod  will  not  only  eliminate  boresight  errors  but  also  free  up  the 
weapons  station  often  occupied  by  the  laser  spot  tracker,  resulting  in  lethality  and 
survivability  benefits. 
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II.  QUANTUM  WELL  STRUCTURE  AND  INTERSUBBAND 

ABSORPTION 


A.  INTRODUCTION 

Intersubband  absorption  is  an  important  aspect  in  the  design  of  an  infrared 
detector  and  has  been  studied  extensively  using  multi-band  methods  (Shik,  1992,  Peng  et 
al.,  1995,  Flatte  et  al,  1995)  and  also  using  band-mixing  (Peng  et  al,  1995).  These 
methods  for  studying  absorption  are  very  sophisticated  and  may  not  be  suitable  for  use  in 
designing  QWIPs.  When  quantum  well  structures  are  formed  using  large  bandgap 
materials,  the  envelope  function  approximation  can  be  used  instead  to  provide  a  more 
simple  approach  to  determining  design  parameters.  In  the  envelope  function 
approximation,  the  Schrodinger  equation  is  solved  by  using  the  transfer  matrix  method, 
which  has  been  discussed  by  several  researchers  (Vassel  et  al,  1983,  Lui  et  al,  1986, 
Wang  et  al,  1989)  and  has  become  the  method  most  commonly  used  to  model  quantum 
well  technology  today.  This  chapter  will  cover  the  transfer  matrix  method,  band 
structure,  selection  rules  for  energy  transitions  and  intersubband  absorption. 


B,  TRANSFER  MATRIX  METHOD 


The  Schrodinger  equation  inside  an  infinite  square  potential  well,  which  is  shown 
in  Figure  2.1,  is  given  by: 


2m^  dz^ 


(2.1) 


where  is  the  mass  of  the  electron,  fi  is  the  reduced  Planck’s  constant,  z  is  the  growth 
direction  and  ij/^  andE’^  are  the  wavefunction  and  energy  eigenvalues,  respectively. 

Finding  the  solution  to  (2.1)  is  relatively  easy  due  to  the  conditions  that  the  wave 
functions  must  go  to  zero  at  the  edges  of  the  well  and  that  the  electron  must  be  within  the 
well.  In  this  case  the  quantized  energy  levels  have  the  values  (West  et  al,  1985): 
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(2.2) 


^n=- — TT  forn=  1,2,3... 

2m/ 

where  L  is  the  width  of  the  well.  The  wavefunetions  in  this  ease  are  given  by 


[2  . 

f  TTnz^ 

—  sm 

u 

1  L  J 

(2.3) 


where  the  ,  |—  term  is  known  as  the  normalization  faetor  and  eomes  from  the  eondition 
L 

that 


j  \¥ni^)  f  dz  =  \. 


(2.4) 


Figure  2.1  Sehematie  drawing  of  an  infinite  square  quantum  well  with  energy  levels. 


In  the  ease  of  an  infinite  potential  well  the  approaeh  used  above  is  relatively 
simple,  but  it  beeomes  tedious  and  diffieult  to  solve  when  multiple  layers  are  added  as  in 
the  ease  of  a  finite  step  quantum  well  (Yuh  et  al,  1989).  In  this  seetion  the  transfer 
matrix  method  will  be  diseussed  and  shown  to  be  able  to  handle  potential  wells  with 
many  layers. 

In  a  heterostrueture,  the  quantum  well  is  formed  by  eombining  semieonduetor 
materials  that  have  differing  bandgaps.  The  simplest  form  of  this  is  a  square  quantum 
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well  in  which  material  A  is  surrounded  by  material  B,  which  has  a  higher  bandgap  and 
which  is  thick  enough  not  to  allow  for  tunneling  of  the  confined  electron.  In  this  case  the 
electrons  will  be  confined  to  the  smaller  bandgap  region  and  using  the  envelope  function 
approximation  (also  known  as  the  effective  mass  approximation)  the  energy  eigenvalues 
can  be  derived  (Kane,  1957). 

The  envelope  function  approximation  makes  use  of  two  assumptions;  (1)  that  the 
potential  difference  between  material  A  and  B  at  the  interface  is  localized  (i.e.  step 
function)  and  (2)  the  interface  potential  does  not  mix  the  band-edge  wavefunctions,  but 
shifts  them.  With  these  assumptions  the  wavefunctions  take  the  form  (Coon  et  al,  1984) 

T(r)  =  M(r)e''^-'’<"’V(2)  (2.5) 


where  ^(r)  is  the  bandedge  cell  periodic  function,  z  is  the  crystal  growth  direction, 

^ikp(x,y)  -g  plane-wave  perpendicular  to  z,  p(x,y)  is  the  direction  vector  in  the  xy 
plane,  k  is  the  transverse  wave  vector,  and  y/’{z)  is  the  envelope  wave  function  that 
satisfies  the  Schrodinger  equation 


2m  (z)  dz‘ 


-  +  V(z) 


y/{z)  =  Ey/{z) 


(2.6) 


where  m* {z)  is  the  effective  mass  in  the  direction  of  growth  and  V [z)  is  the  potential 
height. 


The  transfer  matrix  method  utilizes  the  piece-wise  constant  potential  V [z)  to 

simplify  much  of  the  calculations  of  the  envelope  function  approximation  so  that  the 
wave  functions  in  each  layer  are  given  as  a  linear  combination  of  a  right-traveling  wave, 
,  and  a  left-traveling  wave,  ,  where  is  the  wave  vector  in  the  n*  region  (Note 
that  is  real  in  the  layer  where  E^>V  (z)  and  imaginary  otherwise).  Using  this  fact 
and  applying  the  boundary  conditions 


¥ni^n)  =  ¥nA^n) 


1  ^¥ni^n) 


1 


(2.7) 
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we  can  mathematically  determine  the  energy  eigenvalues. 


V(z) 


Figure  2.2  Arbitrary  potential  well  of  a  multi-layered  quantum  well  structure. 


Figure  2.2  shows  an  arbitrary  quantum  well  structure  in  which  the  thickness  of 
each  layer  is  given  by  and  m*  is  the  effective  mass  of  the  electron  in  each 

region.  The  wavefunctions  for  each  region  are  given  by 


-ik„z 


n  =  0,l,2,...,n,n+l 


(2.8) 


where  and  are  the  coefficients  of  the  left  and  right  traveling  wave  in  the  n*’’  region 
and  the  wave  vector  in  the  n'*'  region  is  defined  as 


(2.9) 


where  is  the  potential  height  in  the  n**'  region.  By  applying  the  boundary  conditions  a 

relationship  can  be  developed  between  the  coefficients  of  the  wavefunctions  in 
neighboring  regions. 

n=  1,2,...,N-1  (2.10) 


^n+l 

=  M„ 

n 

where  is  the  transfer  matrix  for  the  interface  at  z^ 
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(2.11) 


M 


n 


1  (1+X,)e'‘*" 


(l+Z,)e“'^*"“*"+‘^^" 


y  *.'"..L 

/C«  ,  * 


(2.12) 


From  this  formula  a  relationship  can  further  be  found  between  the  first  and  the  last  region 
in  the  quantum  structure 


A 


=  M^_j...M2Mi 


=  M 

mil 

m,2 

~A 

_m2i 

m22_ 

A. 

(2.13) 


where  M  is  the  composite  transfer  matrix  for  the  entire  structure  and  is  given  by 


M  =  M^_i...M2Mi 


m,i  mi2 


(2.14) 


If  the  potential  heights  in  the  two  outer  regions  are  higher  than  those  in  between  them 
there  will  be  bound  energy  states.  For  these  bound  states  the  wavefunctions  must  go  to 
zero  at  ±oo  and  this  requirement  means  that  =0  and  Bjy  =  0  : 


1 

1 

1 

s 

mi2 

1 

o 

1 _ 

1 

O 

1  _ 

_ J 

- 1 

<N 

1 - 

> 

1 _ 

which  is  satisfied  when 


(2.15) 


m22(B)  =  0  (2.16) 

Equation  (2.16)  is  then  tantamount  to  an  eigenvalue  condition  on  the  allowed  energies. 
We  wrote  a  program  to  plot  m22{E)  versus  E  and  locate  the  zeroes  of  the  function,  which 
determine  the  energy  eigenvalues.  The  zeroes  of  m22{E)  are  found  by  searching  for 
values  of  E  such  that  m22{E)\m22{E  +  EE)\  <  0  ;  in  this  way  we  know  that  the  root  lies 

between  E  and  E  +  EE  .  The  accuracy  of  the  eigenvalue  is  determined  by  EE ;  for  most 
cases  a  value  of  EE  =  1  meV  was  used. 
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An  example  of  the  transfer  matrix  method  for  a  square  quantum  well  is  shown  in 
the  next  seetion,  where  the  energy  levels  for  the  quantum  well  shown  in  Figure  2.3  are 
derived. 

V(z) 


m 


V. 


m 


V2  =  0 


V3=V, 


m*3=m*^ 


kg  k.| 


z 


0 


Figure  2.3  Sehematie  diagram  of  a  finite  square  well  with  parameters  used  for 

ealculation. 


The  eomposite  transfer  matrix  M  for  this  quantum  well  ean  be  written  as 


M  =  M2M,=i 


l+Zi  l-Z 
I-Zl  I+Zl 


(2.17) 


2m  2m 

where  k^=k^=  iK  =  i.  — ~  k2  =  ^  ,  so  that, 

\  h  \  h 


m 


22' 


i{k2+k2)a 


^-i{k2-k2)a 


(2.18) 


The  bound  states  oeeur  when  m22iE)=0.  Figure  2.4  shows  a  plot  of  ^  funetion 

of  energy. 


12 


m22(E)  vs.  Energy 


Figure  2.4  Plot  of  versus  E  for  InP ! In^^^Ga^^As  well  with  a  =  150  A. 

Although  only  the  conduction  band  was  considered  in  this  case,  the  valence  band  follows 
the  same  procedure  and  the  energy  eigenvalues  in  the  valence  band  can  be  solved  for 
with  the  same  program  keeping  in  mind  that  the  effective  masses  and  potential  heights 
will  differ  in  the  valence  band.  Figure  2.5  shows  the  conduction  band  energy  levels  and 
wavefunctions  of  the  quantum  well  in  Figure  2.4. 
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Conduction  Band  Energy  Levels  and  Wavefunctions 


Figure  2.5  Conduction  band  energy  levels  and  wavefunctions  of  the  quantum  well  in 

Figure  2.4. 

C.  SELECTION  RULES  FOR  OPTICAL  TRANSITIONS 

In  a  quantum  well  structure  made  of  semiconductors,  the  wavefunctions  consist  of 
envelope  and  Bloch  functions.  This  means  that  transitions  between  both  Bloch  and 
envelope  states  are  allowed  within  the  structure.  Transitions  between  Bloch  states 
correspond  to  interband  transitions,  while  transitions  between  envelope  states  correspond 
to  intersubband  transitions.  The  transition  rate  for  either  interband  or  intersubband 
transitions  is  given  as  (Coon  et  al,  1984) 

W  =  — Xl<  I  I >1'  S[E^-E,-no)\  (2.19) 

n  ij 

where  E.  andE^  are  the  initial  and  final  energy  states  and  hO)  is  the  incoming  photon 

energy.  The  initial  and  final  wavefunctions  in  the  transition  rate  are  of  the  form  of 
Equation  (2.5), 
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(2.20) 


The  interaction  potential  is  defined  in  the  dipole  approximation  as  (Coon  et  al,  1984) 


in 


J/2 


Is^n^CiX: 


A. 

e  p 


(2.21) 


where  /  is  the  incident  photon  flux,  c  is  the  speed  of  light,  e  is  the  electron  charge,  £q  is 
the  electric  permittivity,  is  the  index  of  refraction,  a>  is  the  angular  frequency  of  the 

incoming  photon,  e  is  the  unit  vector  of  the  electric  polarization  of  the  incident  light  and 
p  is  the  momentum  operator.  The  integral  in  the  transition  rate  is  known  as  the  transition 
matrix  element,  which  is  defined  as 

=< I I  >=<  t//r)(/4/z) I I  > 

=<  M^(r)e I  Vp  I  ><  y/f{z)  \  Y^iz)  >  (2.22) 

+  <  I  Vp  I  >< M^(r)  | M,(r)  > . 


For  the  case  of  interband  transitions  ( y/^  =  ,  W f-¥c)  second  term  in  the  transition 

matrix  element  is  zero  due  to  the  orthogonality  of  the  Bloch  functions  in  the  conduction 
and  valence  bands.  On  the  other  hand,  the  first  term  is  zero  for  intersubband  transitions 
( \|/^  =y/y,  yf f  =  yf^)  due  to  the  fact  that  the  Hamiltonian  used  to  obtain  y/^  and  y/^  is 

Hermitian.  Taking  this  into  account,  the  transition  matrix  elements  for  both  transitions 
are  as  follows 


M,,  =<  M,(r)  I  Vp  I  M.(r)  X  y/^z)  \  y/^z)  > 
Mi2  =<  1^2 (^)  I  Vp  I  y/,{z)  ><  u^{y)  I  Mi(r)  > 


(2.23) 


where  the  term,  giving  a  delta  function  for  the  conservation  of  the  xy  plane 

momentum,  has  been  dropped.  It  is  easy  to  see  that  the  only  way  to  obtain  non-zero 
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transition  matrix  elements  is  for  both  <  y/^{z)  \  0  and  <  y/y{z)  \  \  0  • 

For  a  square  quantum  well  this  implies  that  only  even-to-even  or  odd-to-odd  interband 
transitions  are  possible,  while  for  intersubband  transitions,  oc  and  only  even-to- 
odd  or  odd-to-even  transitions  are  possible  as  shown  in  figure  2.6. 


Figure  2.6  Sehematic  diagram  of  a  square  quantum  well  showing  allowed  and 

forbidden  (X)  transitions. 

If  a  step  is  added  within  the  well  in  Figure  2.6  then  the  wavefunctions  are  no 
longer  only  even  or  odd,  but  a  combination  of  the  two.  In  this  case  the  transition  matrix 
elements  will  be  non-zero  for  all  transitions  as  is  shown  in  figure  2.7. 
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Figure  2.7 
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Schematic  diagram  of  a  step  quantum  well  showing  that  all  transitions  are 

allowed. 


This  step  is  vital  to  this  thesis  because  it  allows  for  the  excitation  of  electrons  from 
the  and  states  to  the  E^  state  to  be  collected  as  photocurrent.  The  optimization  of 
transition  wavelengths  and  the  transition  probability  are  discussed  in  the  next  section. 


D.  INTERSUBBAND  ABSORPTION 

Absorption  is  an  important  factor  in  the  design  of  a  multicolor  infrared  detector 
because  it  is  one  of  the  only  measurable  quantities  of  a  sample.  The  absorption 
coefficient  is  proportional  to  the  amount  of  incident  light  that  is  absorbed  in  the  sample 
and  therefore  should  be  high  in  order  to  achieve  an  efficient  detector.  This  thesis  will 
include  absorption  from  both  bound-to-bound  and  bound-to-continuum  transitions,  which 
are  shown  in  figure  2.8. 


17 


ho) 

_ A _ / 

_ A _ 

- ! 

hco 

^ 

r - 

(a)  (b) 

Figure  2.8  Diagram  of  (a)  bound-to-continuum  and  (b)  bound-to-bound  transitions. 


For  the  case  of  bound- to-bound  transitions  there  is  only  one  initial  state  and  one 
final  state,  but  there  are  N  =  N^LA  electrons  occupied  in  the  initial  state  due  to  freedom 
of  motion  of  the  electrons  in  the  xy  plane,  where  N  is  the  number  of  electrons  in  the 
initial  state,  is  the  doping  concentration  in  electrons  per  volume,  L  is  the  length  of 

the  well  and  A  is  the  area  of  the  well.  Since  all  of  these  electrons  have  the  same  energy 
along  the  growth  direction  the  sum  in  Equation  (2.19)  is  over  just  one  initial  and  final 
state  multiplied  by  N  .  The  delta  function  in  Equation  (2.19)  can  be  approximated  due  to 
finite  lifetime  of  the  excited  state  by  the  normalized  Eorentzian  function  given  as  (Yariv, 
1997,  pp.  162) 


- (52) - ^ 


(2.24) 


where  E  is  the  full  width  half  maximum  value  of  the  electron  in  the  state  .  Elsing  the 
interaction  potential  from  Equation  (2.21)  and  the  momentum  operator 


(2.25) 


the  transition  rate  can  be  found  as 


'V,-,  =^(IAL) \< 1 4- 1  >1’  S(*®) 


(2.26) 
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where  <  |  |  'P,  >~<  |  |  >  for  intersubband  transitions  and 


C 


B-B 


N,e"n^ 

[m*)  e^n^c 


(2.27) 


For  the  ease  of  bound-to-eontinuum  transitions  there  are  still  N  initial  states  with 
energy  E. ,  but  now  there  is  also  a  eontinuum  of  final  energy  states  and  the  summation 
over  these  final  states  amounts  to  replaeing  the  delta  funetion  with  the  joint  density  of 
states  p{E^),  whieh  is  given  by  (Choi,  1997,  pp.  131-133) 


27r^h\E.-V,) 


(2.28) 


where  ml  is  the  barrier  effeetive  mass  and  FJ,  is  the  barrier  potential  height.  Therefore 
the  transition  rate  for  bound-to-continuum  states  is  defined  as 


W, 


B-C 


2m, 


ihco)\E.-V,) 


\<  ¥f 


A. 

dz 


¥i  >r 


(2.29) 


where 


2+,2 


/  =*=  \2 

2(m  j  epi^c 


(2.30) 


The  absorption  eoefficient  a{fiCO)  is  defined  as 


a{hQ))  = 


number  of  transitions  per  unit  volume  and  time 
ineident  photon  flux 


(2.31) 


and  the  absorption  eoefficient  for  both  bound-to-bound  and  bound-to-continuum 
transitions  can  be  derived. 


cc^.sinco)  =  \<¥f\^\  ¥i  >f  gi^O)) 

nO)  dz 


2m 

a„  AfiO))  =  C„  r  I - ^ — - \<¥f\  —  I ¥i  >1 


d  I  i2 

—  I V/,  >1 

dz 


(2.32) 
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An  example  of  the  absorption  coefficient  for  bound-to-bound  and  bound-to-continuum 
transitions  is  shown  in  figures  2.9  and  2.10.  The  sharp  peak  in  figure  2.9  results  from  the 
bound-to-bound  transitions  occurring  whenever  the  incident  photon  has  energy  equal  to 
E^-E^,  whereas  the  shape  of  the  absorption  curve  in  figure  2.10  is  relatively  broad  due 

to  the  continuum  of  energies  above  the  well.  The  amplitude  of  the  bound-to-continuum 
absorption  coefficient  decreases  as  the  photon  energy  increases  because  the  density  of 
states  decreases  the  further  above  the  well  you  go. 


Bound-to-Bound  Absorption 


photon  energy  (eV) 

Figure  2.9  Bound-to-bound  absorption  curve  for  (Ino.53Gao.47As)o.i5(InP)o.85/ 
Ino.53Gao.47As  /(Ino.53Gao.47As)o.45(InP)o.55  with  a=25  A  and  b=44  A  with  no  external  bias 

and  r  =  10  meV  . 
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Bound-to-Continuum  Absorption 


Figure  2.10  Bound-to-continuum  absorption  curve  for  the  quantum  well  in  Figure  2.9 

with  no  external  bias. 


E.  SUMMARY 

In  this  chapter  the  transfer  matrix  method  was  developed  and  used  to  calculate  the 
energy  levels  and  wavefunctions  in  a  quantum  well.  Selection  rules  for  interband  and 
intersubband  transitions  and  their  importance  in  detector  design  were  discussed  and  the 
absorption  coefficient  for  both  bound-to-bound  and  bound-to-continuum  transitions  was 
derived  theoretically  using  Fermi’s  golden  rule  (2.19).  In  the  next  chapter,  we  will 
discuss  some  of  the  design  aspects  that  are  considered  when  designing  a  multicolor 
QWIP  and  how  the  addition  of  an  external  bias  will  effect  these  considerations. 
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III.  DESIGN  OF  A  QW  STRUCTURE  FOR  INFRARED 

DETECTION 


A.  INTRODUCTION 


As  discussed  in  Chapter  II,  surrounding  a  low-bandgap  material  with  a  high- 
bandgap  material  will  form  a  quantum  well.  The  choiee  of  which  materials  to  use  is 
important  as  it  will  affect  the  transition  wavelengths  as  well  as  the  absorption  coefficient, 
and  therefore  the  efficiency  of  the  detector.  Some  of  the  factors  to  consider  in  the  choice 
of  materials  to  use  are  bandgap,  lattiee  eonstant,  and  what  wavelengths  are  to  be  detected. 
During  this  thesis  multiple  materials  were  tried  in  order  to  optimize  the  structure  so  that  it 
would  deteet  a  wavelength  of  1.06  pm  and  a  wavelength  of  8-10  pm.  No  matter  whieh 
materials  were  used  in  the  structure,  their  bandgaps  and  effective  masses  were 
determined  by  Vegard’s  Law  with  a  bowing  parameter  to  take  into  aecount  non¬ 
uniformities  of  alloying  (Singh,  1993,  pp.  185). 


=  xE^  +  (1  -  x)E^  +  x(l  -  x)C 


AB 


m. 


^  X  1-x^ 

ml  j 


(3.1) 


where  x  is  the  molar  pereentage  of  material  A  and  is  the  experimentally  determined 
bowing  parameter.  An  example  of  how  the  eomposite  energy  gaps,  at  the  L -point ,  and 
relative  effective  masses  change  with  molar  pereentage  is  shown  in  Figures  3.1  and  3.2 
where  the  energy  gaps  and  effeetive  masses  for  eaeh  material  are  given  in  Table  3.1 
(Note  that  the  bowing  parameter  is  not  considered  in  these  cases). 
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Energy  Gap  (eV) 


Table  3. 


Material 

Energy  Gap  (eV) 

Relative  effective  mass  [m* 

AlAs 

2.671 

0.15 

GaAs 

1.424 

0.067 

In  As 

0.36 

0.028 

Energy  gap  and  relative  effeetive  mass  values  for  speeific  materials.  From; 


Singh,  1993,  pp. 184-185 
Energy  Gap  vs.  Molar  Percentage 


Figure  3.1  Diagram  of  Energy  Gap  vs.  Molar  percentage  for 
Alfia^_^As  and  Ga.  Jn^_^As  at  the  T-point. 


24 


Relative  Effective  Mass  vs.  Molar  Percentage 


Figure  3.2  Diagram  of  Relative  effeetive  mass  vs.  Molar  pereentage  for 

Alfia^_^As  and  GaJn^_^As 


The  potential  heights  in  the  conduction  and  valence  band  depend  on  the  difference 
in  energy  gaps  between  two  regions  and  are  not  the  same,  as  might  be  expected,  but  differ 
so  that  the  conduction  band  potential  height  is  larger  than  the  valence  band.  This 
difference  is  difficult  to  predict  theoretically  and  can  only  be  determined  for  each 
material  experimentally,  however,  when  modeling  a  quantum  well  structure  it  is  assumed 
that  to  a  good  approximation  the  potential  heights  in  the  conduction  and  valence  bands  is 
given  as 


In  early  QWIP  research,  studies  focused  on  quantum  well  structures  that  utilized 
bound-to-bound  transitions  (Levine  et  al.  1987).  In  a  bound-to-bound  QWIP  the 
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electrons  in  the  lowest  state  are  excited  by  incoming  photons  into  the  excited  state  and 
then  tunnel  through  the  barrier  under  an  external  bias  to  form  a  photocurrent  as  is  shown 
in  figure  3.3. 


Figure  3.3  Schematic  Diagram  showing  (a)  bound-to-continuum  and  (b)  bound-to- 

bound  electron  transport. 

The  barrier  thickness  was  designed  to  be  large  enough  to  prevent  tunneling  from  the 
lowest  state  (typically  about  500  A);  however,  this  barrier  thickness  reduced  the 
photoresponse  due  to  the  reduction  in  tunneling  probability  of  the  photoelectrons. 
Another  problem  with  these  early  designs  was  that  the  photoresponse  spectrum  was 
relatively  narrow  due  to  there  being  only  two  states  and  was  not  practical  for  infrared 
imaging.  These  problems  ushered  in  the  second  generation  of  QWIPs,  which  used 
bound-to-continuum  transitions  (Hasnain  et  al,  1989).  In  these  structures,  electrons  in 
the  ground  state  are  excited  into  the  continuum  where  they  contribute  to  the  photocurrent 
without  tunneling  (Figure  3.3)  and  therefore  the  barrier  thickness  has  less  bearing  on  the 
photocurrent.  Due  to  the  continuum  of  states  and  the  multiple  transitions  that  can  occur, 
the  photoresponse  spectrum  was  much  wider  than  in  the  bound-to-bound  design,  but  the 
absorption  coefficient  was  much  smaller  due  to  the  low  density  of  states  above  the  well. 
It  is  easy  to  see  that  there  is  a  trade-off  between  having  good  transport  of  photoelectrons 
and  strong  absorption  and  that  both  are  important  in  the  design  of  a  QWIP.  In  the 
following  sections,  these  factors  will  be  used  to  motivate  the  design  of  a  quantum  well 
structure  for  two-color  infrared  detection. 
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B,  STEP  QUANTUM  WELL  STRUCTURE  FOR  DUAL  COLOR 

DETECTION 

In  the  past,  attempts  have  been  made  to  deteet  multiple  “eolors”,  or  wavelengths, 
with  an  infrared  imager  by  using  square  quantum  wells  (Koek  et  al,  1992,  Grave  et  al, 
1992,  Liu  et  al,  1993),  but  as  was  shown  in  seetion  2.3,  transition  seleetion  rules  prohibit 
the  transition  from  the  first  valenee  level  to  the  seeond  eonduetion  level  and  therefore  that 
wavelength  eannot  be  deteeted.  This  implies  that  only  one  wavelength  of  light  ean  be 
detected  with  a  square  quantum  well.  One  method  used  to  detect  more  than  one 
wavelength  of  light  using  square  quantum  wells  was  to  stack  multiple  wells  on  top  of 
each  other  and  hold  each  well  responsible  for  detecting  one  specific  wavelength  of  light 
(Koek  et  al,  1992,  Grave  et  al,  1992,  Liu  et  al,  1993).  This  type  of  structure  is  easy  to 
design,  but  very  difficult  to  fabricate  due  to  the  need  for  multiple  contacts  at  each  well. 
Recently  researchers  have  been  exploring  the  use  of  asymmetric,  or  step,  quantum  wells 
(Ting  et  al,  1997)  to  detect  multiple  wavelengths. 

In  a  step  quantum  well,  the  symmetry  of  the  wavefunctions  is  broken  and 
transitions  that  were  forbidden  in  a  square  quantum  well  can  occur.  For  this  thesis,  the 
optimized  design  of  the  step  quantum  well  is  to  have  two  energy  levels  in  the  conduction 
band,  one  above  and  below  the  step,  which  account  for  the  middle  infrared  wavelength, 
and  one  energy  level  in  the  valence  band,  below  the  step,  that  accounts  for  the  1 .06  pm 
transition  as  illustrated  in  figure  3.4.  The  second  energy  level  in  the  conduction  band  is 
designed  to  be  as  near  the  top  of  the  quantum  well  as  possible  in  order  to  get  a  good 
photocurrent  with  very  little  bias. 
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V(z) 


The  application  of  the  bias  will  cause  the  quantum  well  to  “slant”  and  will  force 
the  second  energy  level  into  the  continuum  allowing  the  photoexcited  electrons  to  be 
extracted  without  tunneling  through  the  barrier.  The  energy  levels  will  be  affected  by  the 
bias  and  therefore  the  peak  detection  wavelengths  will  shift.  The  asymmetry  of  the 

cFb 

quantum  well  will  give  rise  to  a  linear  Stark  shift  which  is  approximately  given  by 

where  e  is  the  electron  charge,  F  is  the  applied  electric  field  in  volts  per  meter  and  b  is 
the  length  of  the  step  in  the  quantum  well.  The  accurate  shift  in  detection  wavelength  is 
important  in  the  design  since  it  is  necessary  to  apply  an  external  bias  to  a  device  in  order 
to  extract  the  photocurrent.  Due  to  the  necessity  of  using  an  external  bias  it  becomes 
important  to  be  able  to  model  how  the  peak  detection  wavelengths  will  change  in  order  to 
optimize  the  QWIP  structure,  in  particular  for  the  1 .06  |im  wavelength.  The  next  section 
will  introduce  a  transfer  matrix  method  for  a  quantum  well  under  a  bias  and  this  method 
will  be  used  to  model  the  change  in  transition  wavelengths  due  to  an  external  bias. 
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C.  TRANSFER  MATRIX  METHOD  UNDER  AN  APPLIED  BIAS 

When  a  quantum  well  structure  is  subjected  to  an  applied  bias  the  potential  well 
will  tilt  in  the  direction  opposite  the  electric  field  and  therefore  the  potential  of  an 
electron  in  the  quantum  well  is  given  as 

V{z)  =  V„-eFz  (3.3) 


where  is  the  potential  height  of  the  well,  eis  the  electron  charge,  F  is  the 

applied  electric  field  in  volts  per  meter  and  z  is  the  direction  of  growth.  An  example  of 
the  effect  of  an  external  bias  on  a  square  quantum  well  structure  is  shown  in  Figure  3.5 


Figure  3.5  Diagram  of  a  square  quantum  well  under  an  external  bias. 
Substituting  Equation  (3.3)  into  the  Schrodinger  equation 


^  2m*  dz^ 


+  V(z) 


(3.4) 


the  following  formula  can  be  derived 


dz' 


where 


(3.5) 
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(3.6) 


Vn  = 


E-Vn 

e 


and  where  is  the  potential  height  in  the  n*  region.  Equation  (3.5)  ean  further  be 
simplified  using  the  following  eoordinate  transformation 

=  (3.7) 


and  the  result  is  of  the  form  of  the  Airy  differential  equation 


(3.8) 


with 


f;j2^2V3 


(3.9) 


The  solution  to  Equation  (3.8)  is  a  linear  eombination  of  Airy  functions  given  as 
(Vatannia  et  al,  1996) 

(z)  =  AJi{p^  (z))  +  B„Bi{Pn  (^))  (3.10) 


where  Ai  and  Bi  are  Airy  functions  of  the  first  and  second  kind,  respectively.  By 
substituting  Equation  (3.10)  into  the  boundary  conditions  given  by  Equation  (2.7)  a 
relationship  can  be  developed  between  the  coefficients  of  the  wavefunctions  in  each 
region. 


=  M„ 

A+i 

Bn'_ 

_Bn+\  _ 

for  n  =  1,2,3,...A-1 


(3.11) 


Eurthermore,  a  relationship  can  be  developed  between  the  coefficients  in  the  first  and  last 
region. 


A 


=  M 

_ 

mil 

m,2 

Bn} 

Bn} 

^21  _ 

Bn 

(3.12) 
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where 


M„  =;r 


Aiia„JBi'ia„)-  a„Ai{a„^^)Bi{a„)  Biia„JBi'ia„)-  CT„Bi{a„)BiXa„j' 
yanAi(an)AiXanJ - AiXaJAiia^J  aJi{aJBiXa^J -  AiXaJBi{a„^^) 


(3.13) 


%+i^y 


and 


C(„ 


V  y 


'iml.eF''''’ 


(3.14) 


/  *  \ 


(3.15) 


The  primes  in  Equation  (3.13)  represent  the  derivatives  of  the  Airy  functions  and  the 

factor  of  ^  comes  from  the  fact  that  [Ai{aJBiX(X„)-  AiXo(JBi(aj)  =  — .  In  the  case  of 

n 

bound  states,  the  wavefunctions  must  decay  to  zero  as  z  ^  ±oo  .  This  constraint  forces 
=  0  since  the  second  order  Airy  function,  Bi(aJ  ,  tends  to  oo  as  z  — >  -oo  .  Similarly, 

Aj^  =0  since  Ai{a^)  tends  to  0°  as  z  — >  0° .  These  conditions  cause  Equation  (3.12)  to 
take  the  form 


1 

1 

1 

_ 1 

1 

0 

1 _ 

0 

1  _ 

<N 

_ 1 

m22_ 

1 - 

to 

1 _ 

(3.16) 


and,  just  like  the  derivation  of  the  unbiased  transfer  matrix  method,  the  requirement  for 
eigen  states  is 

m^2iE)  =  0.  (3.17) 


We  wrote  a  program  to  find  and  plot  m22{E)  versus  E  and  find  the  energy  levels  when 

Equation  (3.17)  is  satisfied.  Eigure  3.6  shows  m22{E)  as  a  function  of  energy  and  figure 

3.7  shows  the  calculated  energy  levels  and  wavefunctions  for  the  conduction  band  of  a 
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step  quantum  well  under  an  external  bias.  The  parameters  of  the  quantum  well  were 
taken  from  a  previously  published  work  (Karunasiri  et  ai,  1990)  so  that  the  results  found 
using  this  program  eould  be  eompared  with  experimental  data. 


m22(E)  vs.  Energy 
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Figure  3.6 


Plot  of  m22(E)  vs.  E  for  Aloji4GaQ^f,As/GaAs/ Alo  isGaog2As  well  with  a=60  A 
and  b=90  A  and  an  external  eleetrie  field,  =  3.5E///m  . 
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Conduction  Band  Energy  Levels  and  Wavefunctions 

T - 1 - 1 - 1 - 1 - 1 - 1 - r 
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z  (angstroms) 


Figure  3.7  The  first  three  conduction  band  energy  levels  and  wavefunctions  of  the 

quantum  well  in  figure  3.6. 

The  eigen  states  and  wavefunctions  of  the  step  well  formed  in  the  valence  band 
can  be  obtained  using  the  same  approach  as  above  with  corresponding  potential  heights 
and  effective  masses.  Figure  3.8  shows  the  calculated  energy  levels  and  wavefunctions 
in  the  valence  band  for  heavy  holes. 
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Valance  Band  Energy  Levels  and  Wavefunctions 


Figure  3.8  Valence  band  energies  and  wavefunctions  for  the  step  well  in  Figure  3.6. 

The  introduction  of  the  external  bias  will  cause  a  shift  in  the  energy  levels  of  the 
asymmetric  quantum  well  and  therefore,  the  absorption  coefficient  defined  by  Equation 
(2.32)  will  also  change.  For  the  QWIP  design  used  in  this  thesis,  the  addition  of  the  bias 
will  force  the  excited  state,  which  was  very  near  the  top  of  the  well,  into  the  continuum 
so  that  only  the  bound-to-continuum  absorption  coefficient  is  considered  when  a  bias  is 
introduced.  Figure  3.9  demonstrates  how  the  external  bias  will  affect  the  peak  detection 
wavelength. 
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Bound-to-Continuum  Absorption 


wavelength  (|j.m) 

Figure  3.9  Plot  of  the  bound-to-eontinuum  absorption  vs.  wavelength  for  the 
(^'*0.53^^0.47^'^)o.l5(^'^P)o.85  ^  KI>^0.53^‘^OAl^^)o.55i^’^^)oA5  Well  wlth  a  —  25  A  and  b  — 

44  A  and  with  a  bias  of  0  V/|im  and  1.5  V/|im. 

This  section  has  introduced  a  transfer  matrix  method  for  a  quantum  well  structure 
under  an  external  bias.  This  method  was  used  to  calculate  the  energy  levels, 
wavefunctions  and  absorption  coefficient  for  a  quantum  well,  which  is  useful  in 
optimizing  detector  design.  The  next  section  will  compare  the  output  of  the  computer 
model  designed  in  this  thesis  with  experimental  data  in  order  to  tell  the  accuracy  of  the 
computer  model. 


D.  COMPARISON  WITH  EXPERIMENTAL  DATA 

The  quantum  well  structure  that  the  computer  model  was  compared  to  came  from 
a  paper  written  by  Karunasiri  et  al.  (1990).  This  paper  utilized  a 
/  GaAs  I  AIq  ^^GaQ^2^s  step  quantum  well  with  a  =  60  A  and  b  =  90  A.  This 

structure  has  four  energy  levels  in  its  conduction  band,  but  the  paper  focused  on  the 
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bound-to-bound  absorption  from  the  1-2  and  1-3  transitions  due  to  the  minute  absorption 
strength  of  the  1-4  transition.  Figure  3.10  shows  the  energy  levels  of  the  conduction  band 
and  Table  3.2  shows  the  effective  masses  and  energy  gaps  used  in  the  computer  model. 


Energy  Levels  of  the  Asymmetric  Quantum  Well 
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Figure  3.10  Energy  levels  of  the  !  GaAs !  well  with  a  =  60A 

and  b  =  90A. 


Layer 

m*  (kg) 

E,  (eV) 

Alo,44Gao,56As 

0.1082 

1.9727 

GaAs 

0.067 

1.424 

Alo,i8Gao,82As 

0.0794 

1.6485 

Table  3.2  Parameters  used  in  the  quantum  well  in  Figure  3.10 

Since  this  is  a  step  quantum  well  all  energy  transitions  are  allowed  and  only  the  bound- 
to-bound  absorption  is  considered  in  this  analysis  to  compare  with  the  experimental  data. 
Figure  3.11  shows  the  calculated  absorption  coefficient  for  the  1-2  and  2-3  energy 
transitions  under  both  positive  and  negative  external  bias. 
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Bound-to-Bound  Absorption  for  1-2  and  1-3  Transitions 


Figure  3.11  Bound-to-bound  absorption  curves  for  1-2  and  1-3  transitions  of  the 
quantum  well  in  Figure  4.6  with  r=15meF  . 

Since  the  experimental  data  is  measured  in  absorbance,  the  absorption  coefficient  is 
converted  to  absorbance  using  the  following  relation 

Absorbance  =  x{\oge)a{ho))  =  W  (\l°)){\oge)a{hcd) 

(3.18) 

=  (2.56*10"'  cm)(log  e)a{n  a) 

where  W  is  the  width  of  the  absorbing  material  and  the  sin^(17°)  term  comes  from  the 

polarization  of  the  light  hitting  the  device  at  the  Brewster  angle  {6  =  13°).  Using  this 
equation  the  absorbance  can  be  found  for  both  transitions  and  is  shown  in  Figure  3.12. 
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Bound-to-Bound  Absorbance  for  1-2  and  1-3  Transitions 


Figure  3.12  Bound-to-bound  absorbance  curves  for  1-2  and  1-3  transitions  of  the 
quantum  well  in  Figure  4.6  with  r=15meF  . 


This  Figure  can  be  compared  with  the  absorbance  curves  obtained  experimentally,  which 
is  shown  in  Figure  3.13  (Note  that  the  absorbance  data  in  the  Figure  is  plotted  in 
decreasing  order  of  wavenumber). 
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Figure  3.13  Bound-to-bound  absorbance  curves  obtained  experimentally  using  the 
quantum  well  in  Figure  4.6.  From:  Karunasiri  et  a/.,  1990. 


Comparing  Figures  3.12  and  3.13  it  can  easily  be  seen  that  the  values  for  the  absorbance 
obtained  from  our  model  and  those  obtained  experimentally  are  very  close  and  that  the 
computer  model  used  in  this  thesis  is  a  good  approximation  to  the  real  world. 


E.  SUMMARY 

This  chapter  introduced  some  of  the  parameters  that  affect  the  design  of  a  QWIP 
and  a  method  for  detecting  multiple  infrared  wavelengths  using  an  asymmetric  quantum 
well  design.  The  optimum  design  characteristics  used  in  this  thesis  were  also  presented 
along  with  a  motivation  for  modeling  the  effects  of  an  external  bias  on  a  quantum  well 
structure.  A  transfer  matrix  method  utilizing  Airy  functions  was  derived  to  model  the 
effects  of  an  external  bias  and  this  method  was  used  to  show  how  the  peak  detection 
wavelength  shifts  due  to  the  bias.  The  data  obtained  using  our  computer  model  was  then 
compared  to  experimental  data  and  shown  to  be  a  good  approximation  to  real  life, 
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granting  greater  eonfidenee  in  any  information  obtained  using  the  model.  The  next 
chapter  introduces  the  optimum  design  of  the  QWIP  structure  used  in  this  thesis  along 
with  an  analysis  of  how  it  will  react  to  an  external  bias. 
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IV.  DESIGN  PARAMETERS  OF  THE  TWO-COLOR  DETECTOR 


A.  INTRODUCTION 

As  was  mentioned  in  chapter  III,  many  different  material  systems  were  considered 
in  designing  the  asymmetric  well  used  in  this  thesis.  At  first,  a 
Alfia^_^As  /  GaAs  /  AlyGa^_yAs  step  well  was  considered  due  to  the  abundance  of 

experimental  data  on  these  two  materials  and  because  they  are  lattice  matched,  which 
would  make  them  much  easier  to  grow.  This  structure  could  be  designed  to  detect  the 
middle  infrared  wavelength  (MIRW),  but  could  not  detect  the  1.06  |im  wavelength  due  to 
its  large  bandgap.  This  prompted  the  use  of  materials  with  smaller  bandgaps,  which 
could  be  identified  by  using  Figure  4.1  (Tu  et  al,  1992,  pp.  158). 

Lattice  mismatch,  f 


0.55  0.60  0.65 

Lattice  constant  (nm) 

Figure  4.1  Diagram  of  bandgap  vs.  Lattice  constant  for  different  materials. 
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The  next  -based  quantum  well  design  that  was  tried  eonsisted  of 

GaAs  /  In^Ga^_^As  /  InyGa^_yAs  and  eould  be  designed  to  deteet  both  the  MIRW  and  1.06 

|im  wavelengths,  but  not  at  the  same  time.  The  reason  that  both  wavelengths  eould  not 
be  deteeted  at  the  same  time  using  these  materials  was  that  in  order  to  deteet  the  1.06  |im 
wavelength  the  energy  gaps  had  to  be  widened,  but  this  eaused  the  two  conduetion  band 
energy  levels  to  beeome  closely  spaced,  which  pushed  the  MIRW  outside  of  the  8-10  |im 
design  window.  Another  problem  with  the  design  of  this  quantum  well  was  that  if  the 
molar  percentage  of  Indium  became  large  there  was  a  significant  lattice  mismatch 
between  the  two  materials,  which  would  introduce  dislocations  into  the  sample.  In 
order  to  lessen  the  lattice  mismatch,  a  quantum  well  consisting  of 
InP  !  In^Ga^_^As  /  InyGa^_yAs  was  designed  because  the  energy  gap  of  InP  and  GaAs  are 

about  the  same  and  because  the  two  materials  have  much  closer  lattice  constants  to 
provide  for  easier  growth.  This  quantum  well  structure  could  also  be  designed  to  detect 
the  MIRW  and  1.06  |im  wavelengths,  but  the  problem  remained  that  both  wavelengths 
could  not  be  detected  at  the  same  time  because  the  energy  gap  of  InP  was  still  too  large. 
Following  Figure  4.1  this  left  only  one  choice  of  possible  materials  for  the  final  design: 
the  quaternary  InGaAsP  and  As  .  These  materials  provided  the  advantage  of 

being  lattice-matched  and  also  by  varying  the  amount  of  InP  in  the  quaternary  the 
bandgap  could  be  brought  down  to  the  point  where  both  wavelengths  could  be  detected. 
The  following  section  will  discuss  the  quantum  well  that  was  designed  using  these 
materials  and  will  show  the  theoretical  results  for  the  absorption  with  and  without  an 
external  bias. 


B,  OPTIMUM  QUANTUM  WELL  DESIGN 

One  of  the  problems  with  using  the  quaternary  material  in  the  design  of  the 
quantum  well  was  that  very  little  experimental  data  was  readily  available  and  no 
theoretical  data  could  be  found.  However,  a  paper  written  by  Vurgaftman  et  al.  (2001) 
provided  all  of  the  necessary  experimental  data  not  only  for  the  quaternary  material,  but 
also  for  every  other  material  that  was  used  throughout  this  thesis.  Using  the  values  from 
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this  paper,  an  optimized  quantum  well  was  designed  with  the  following  layer  structure: 
(/Hq  53^6(0 47^45)0  j5(/n/’)og5 //«(,  53^00  47^45 /(/nQ53Gao 47^45)0  55 (/nP),,  45  with  a  thickness  of 
the  well  of  25  A  and  of  the  step  of  44  A.  Figure  4.2  shows  the  optimized  quantum  well 
structure,  energy  levels,  and  transition  wavelengths  and  Table  4.1  shows  the  values  of  the 
effective  masses  and  energy  gaps  in  each  layer  of  the  optimum  design. 


Energy  Levels  of  the  Optimized  Quantum  Well 


0.4 

"  1  1  1  1  1  1  1  1 

0.2 

j 

1  i 

. 

0 

E.,=  0.1576  eV 

E2=  0.2995  eV 

-0.2 

E^.,= -0.8698  eV 

- 

0) 

E^2=  -0.9596  eV 

-G  4 

0) 

E^3=  -1 .0030  eV 

c 

LU 

3,,=  8.739p.m 

-0  6 

^2=  1 .060  p,m 

-0.8 

- 

.1 

-1.2 

_ 1 

_ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 

J _ 

-10  0  10  20  30  40  50  60  70  80 

z  (angstroms) 

Figure  4.2  Optimized  quantum  well  structure  and  energy  levels  with  no  external  bias. 


Layer 

m  (kg) 

Kh  (kg) 

E,  (eV) 

(lno.53Gao.47As)o.i5(lnP)o.85 

0.0676  Wg 

0.4905 

1.3159 

lno.53Gao.47As 

0.0365 

0.3403 

0.8161 

(lno.53Gao.47As)o.55(lnP)o.45 

0.0519 

0.4242 

1.1181 

Table  4.1  Parameters  used  to  design  the  final  quantum  well  structure  used  in  this 
thesis,  mlf^  is  the  heavy-hole  effective  mass  and  is  the  electron  mass. 

It  is  important  to  mention  that  the  values  for  the  energy  gaps  found  in  Table  4.1  are  the 

values  when  the  temperature,  T  =  OK  .  The  reason  for  choosing  these  values  over  the 
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room  temperature  {T  =  290K )  values  is  that  in  order  to  limit  the  dark  current  of  the 
device  under  operating  conditions  it  must  be  cooled  down  io  T  <  11 K  .  The  difference  in 
the  energy  gaps  and  the  energy  levels  of  this  structure  at  7  =  OK  and  T  =  11 K  is  small 
so  that  for  design  purposes  it  was  assumed  that  the  device  is  operating  at  7  =  OK . 
Figures  4.3  and  4.4  show  the  effects  of  external  bias  on  the  quantum  well  structure  and 
energy  levels. 


Energy  Levels  of  the  Optimized  Quantum  Well 


Figure  4.3  Optimized  quantum  well  structure  and  energy  levels  with  Fext=^.svl fim  . 
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Energy  Levels  of  the  Optimized  Quantum  Well 


It  is  easy  to  see  in  these  Figures  that  the  seeond  energy  level  has  indeed  been  pushed  up 
into  the  eontinuum  where  it  ean  easily  ereate  a  photoeurrent. 

The  main  factor  to  consider  in  the  evaluation  of  this  design  is  how  the  application 
of  the  external  bias  will  affect  the  quantum  well’s  ability  to  detect  the  designated 
wavelengths  of  1.06  |im  and  8-10  |im  band.  Figure  4.5  shows  a  comparison  between  the 
absorption  coefficient  for  the  unbiased  well  where  bound-to-bound  absorption  occurs  and 
for  the  case  when  =\.5Vl /ilm  where  bound-to-continuum  absorption  occurs. 
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Bound-to-Bound  and  Bound-to-Continuum  Absorption  vs.  Wavelength 


Figure  4.5  Comparison  of  the  absorption  eoeffieients  for  the  unbiased  well  (bound- 
to-bound)  and  the  biased  well  (bound-to-continuum)  with  Ff,y^i=\5V j ^im  and  r=10ff!eF  . 

This  Figure  demonstrates  the  relative  strengths  of  the  bound- to-bound  and  bound-to- 
eontinuum  absorption  eoeffieients.  The  bound-to-bound  absorption  eoeffieient  has  a 
narrow  speetral  width  eompared  to  the  broad  bound-to-eontinuum  speetrum  as  illustrated 
in  Figure  4.6  using  the  normalized  absorption  eoeffieients. 


46 


Normalized  Absorption  vs.  Wavelength 


wavelength  (|j.m) 


Figure  4.6  Comparison  of  the  normalized  absorption  eoeffieients  for  the  unbiased 
well  (bound-to-bound)  and  the  biased  well  (bound-to-eontinuum)  with  Fgxt=\.5V j /im  and 

r=10weK  . 

Although  the  shift  from  bound-to-bound  to  bound-to-continuum  absorption  has  dropped 
the  relative  height  of  the  absorption  spectrum  it  has  not  dropped  to  a  point  where  it  would 
be  difficult  to  measure.  The  broadening  of  the  absorption  spectrum  has  also  ensured  that 
the  quantum  well  structure  will  still  be  able  to  detect  the  MIRW  reasonably  well  and  as 
long  as  the  external  bias  is  not  too  large,  the  QWIP  should  operate  as  designed.  It  is 
important  to  mention  that  this  thesis  has  not  considered  interband  absorption  due  to  the 
difficulties  in  defining  and  integrating  the  Bloch  functions  in  each  material,  but  it  can 
safely  be  assumed  that  the  interband  absorption  curve  will  follow  the  same  pattern  that 
the  intersubband  absorption  curve  did  and  therefore  the  quantum  well  structure  will  also 
detect  the  1 .06  |im  wavelength  reasonably  well. 
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C.  SUMMARY 

This  chapter  introduced  the  optimized  design  of  the  quantum  well  structure  used 
in  this  thesis  along  with  some  of  the  reasons  why  the  materials  used  in  this  quantum  well 
were  chosen  over  other  materials.  It  was  also  shown  that  the  optimized  quantum  well 
structure  would  operate  as  a  multicolor  detector  for  the  specified  wavelengths  with  or 
without  an  external  bias. 
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V.  CONCLUSION 


An  asymmetric  quantum  well  device  capable  of  achieving  two-color  detection 
through  the  use  of  bandgap  engineering  was  demonstrated  and  the  properties  of  this 
device  were  derived  and  discussed.  The  ability  of  the  optimized  device  to  detect  both 
short  infrared  wavelengths  and  middle  infrared  wavelengths  could  be  useful  as  a  dual 
infrared  imager/laser  spot  tracker  in  military  aircraft  for  the  use  of  accurate  ordinance 
delivery. 

The  detailed  quantum  well  structure  and  eigen  energy  levels  of  the  device  were 
derived  through  the  use  of  a  transfer  matrix  method  with  effective  mass  approximation 
and  the  intersubband  absorption  spectrum  for  both  bound-to-bound  and  bound-to- 
continuum  transitions  was  derived  and  utilized  to  optimize  the  device  parameters.  A 
transfer  matrix  method  using  Airy  functions  was  derived  to  model  the  effects  of  an 
external  bias  on  the  device  and  the  output  of  this  computer  model  was  compared  to 
experimental  data  and  shown  to  be  a  good  approximation  to  real  life.  The  reasoning  for 
the  choice  of  the  materials  used  in  the  optimized  design  was  discussed  including  a 
discussion  of  why  other  materials  that  were  tried  did  not  work  for  the  chosen 
wavelengths. 

This  thesis  has  shown  that  an  efficient  laser  spot  tracker/infrared  imager  quantum 
well  device  is  plausible  and  further  study  is  warranted.  A  sample  of  the  optimized  design 
should  be  grown  and  tested  in  order  to  compare  experimental  results  with  the  theory 
presented  in  this  thesis. 
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APPENDIX.  [PROGRAMMING  NOTES  AND  SELECTED 

PROGRAMS] 


Since  many  of  the  ealculations  presented  in  this  thesis  involved  a  substantial 
amount  of  matrix  operations  the  choiee  of  Matlab  to  create  the  computer  model  was 
manifest.  Matlab  proved  to  be  an  extremely  efficient  and  reliable  tool  for  deriving  the 
energy  levels,  wavefunctions  and  absorption  spectra  for  the  unbiased  quantum  well 
devices.  However,  the  use  of  Matlab  to  model  the  biased  quantum  well  had  some 
limitations  due  to  the  Airy  funetions.  Matlab  contains  the  numerical  values  of  the  Airy 
functions,  but  as  expected,  these  values  become  unmanageable  when  the  argument  in  the 
Airy  functions  becomes  very  large  or  small.  An  attempt  was  made  to  manually  insert  the 
asymptotic  versions  of  the  Airy  functions  into  the  program  to  overcome  this  limitation, 
but  this  too  proved  to  be  unmanageable.  Therefore  the  program  written  to  model  the 
quantum  well  under  an  applied  field  is  limited  in  its  scope.  The  following  is  a  seleetion 
of  programs  used  in  this  thesis  to  obtain  the  theoretical  data  discussed  above. 
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%This  program  is  titled  aqw_InGaAsP% 
%Created:  January  2002 
%Author:  Kevin  Lantz 


%This  is  a  program  to  find  the  energy  levels  in  a  two-step  quantum  well. 

%There  are  three  interfaces  in  this  well  and  therefore  will  be  three  transfer  matrices  to  calculate. 
%The  basic  design  of  the  well  is  as  follows: 

%Region  1  (-Inf->0)->InGaAsP 

%Region  2  (0->a)— >InGaAs 
%Region  3  (a->L)->InGaAsP 
%Region  4  (L->Inf)— >InGaAsP 


clc 

global  E 

%  inputs - 

XI  =  input('Percentage  of  InGaAs  (0-1)  in  region  1  '); 

X2  =  input('Percentage  of  In  (0-1)  in  region  2  '); 

X3  =  input('Percentage  of  InGaAs  (O-I)  in  region  3  '); 

widthi  =  Inf;  %Extends  from  -Infinity  to  0; 
width2  =  input('Length  of  0->a  in  angstroms:  '); 

width3  =  input('Length  of  a->L  in  angstroms:  '); 

width4  =  Inf;  %Extends  from  L  to  Infinity; 


%Energy  gap  values  in  each  region - 

egl  =  1.4236*(1-X1)  -t  0.81612*X1  -  0.13*(X1*(1-X1)); 
eg2  =  1.519*(1-X2)  -t  0.417*X2  -  0.477*(X2*(1-X2)); 
eg3  =  1.4236*(1-X3)  -t  0.81612*X3  -  0.13*(X3*(1-X3)); 
eg4  =  1.4236*(1-X1)  -t  0.81612*X1  -  0.13*(X1*(1-X1)); 

%Potential  Heights  in  Conduction  Band - 

VI  =  0.6*(egl  -  eg2); 

V2  =  0;  %Choose  the  origin  as  the  bottom  left  comer  of  region  2. 

V3  =  0.6*(eg3  -  eg2); 

V4  =  V1; 

% - 

y  =  1 ;  %Initial  value  of  y  to  correct  the  logic  of  the  later  condition; 
a=l; 

%constants - 

hbar  =  1.055e-34;  %hbar  in  ev*m 

mO  =  9.1  le-31;  %Rest  mass  of  the  electron  in  eV 

q=  1.602e-19; 

L 1  =  width  1  *  1  e- 1 0;  %Change  input  of  angstroms->meters 

L2  =  width2*le-10; 

L3  =  (width2-t-width3)*le-10; 

L4  =  width4*le-10; 

%Effective  electron  mass  in  each  region - 

m_effl  =  l/((Xl/(0.0365*m0))  ((l-Xl)/(0.0795*m0))); 

m_eff2  =  l/((X2/(0.026*m0))  ((l-X2)/(0.067*m0))); 

m_efG  =  l/((X3/(0.0365*m0))  ((l-X3)/(0.0795*m0))); 
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m_eff4  =  l/((Xl/(0.0365*m0))  +  ((l-Xl)/(0.0795*m0))); 
for  E  =  0.0001 :0.0001:(V1-0.0001); 


%Determine  the  values  of  k - 

kl  =  sqrt((2*m_effl*(E  -  Vl)*q)/(hbar^2)); 
k2  =  sqrt((2*m_eff2*(E  -  V2)*q)/(hbar^2)); 
k3  =  sqrt((2*m_eff3*(E  -  V3)*q)/(hbar^2)); 
k4  =  sqrt((2*m_eff4*(E  -  V4)*q)/(hbar^2)); 

%Determine  the  transfer  matrices - 


xl  =  0;  %Interface  Locations; 
x2  =  L2; 
x3  =  L3; 

mil  =((1/2)*((1  +((kl*m_eff2)/(k2*m_effl)))*exp(i*(xl*(kl  -  k2))))); 
ml2  =  ((1/2)*((1  -  ((kl*m_eff2)/(k2*m_effl)))*exp(-(i*(xl*(kl  +k2)))))); 
m21  =((1/2)*((1  -((kl*m_eff2)/(k2*m_effl)))*exp(i*(xl*(kl  +k2))))); 
m22  =  ((1/2)*((1  +  ((kl*m_eff2)/(k2*m_effl)))*exp(-(i*(xl*(kl  -  k2)))))); 

m_ll  =  ((1/2)*((1  +  ((k2*m_eff3)/(k3*m_eff2)))*exp(i*(x2*(k2  -  k3))))); 
m_12  =  ((1/2)*((1  -  ((k2*m_eff3)/(k3*m_eff2)))*exp(-(i*(x2*(k2  +  k3)))))); 
m_21  =  ((1/2)*((1  -  ((k2*m_efG)/(k3*m_eff2)))*exp(i*(x2*(k2  +k3))))); 
m_22  =  ((1/2)*((1  +  ((k2*m_efG)/(k3*m_eff2)))*exp(-(i*(x2*(k2  -  k3)))))); 

m_ll_  =  ((1/2)*((1  +  ((k3*m_eff4)/(k4*m_efG)))*exp(i*(x3*(k3  -  k4))))); 
m_12_  =  ((1/2)*((1  -  ((k3*m_eff4)/(k4*m_efG)))*exp(-(i*(x3*(k3  +k4)))))); 
m_21_  =  ((1/2)*((1  -  ((k3*m_eff4)/(k4*m_efB)))*exp(i*(x3*(k3  +k4))))); 
m_22_  =  ((1/2)*((1  +  ((k3*m_eff4)/(k4*m_efB)))*exp(-(i*(x3*(k3  -  k4)))))); 


Ml  =[mll,ml2;m21,m22]; 

M2  =  [m_ll,m_12;m_21,m_22]; 

M3  =  [ml  l_,m_12_;m_21_,m_22_]; 


%Now  we  find  an  expression  for  M22 


B  =  M3*M2*M1; 
J22  =  B(2,2); 


%Program  used  to  calculate  and  plot  the  wavefunctions 


Psiaqw 

%Defme  the  function  y- 


z  =  y; 

y  =  real(J22); 

%Define  condition  for  energy  eigenvalues - 

ifz.*y  <0 
disp(E) 

Eout(a)  =  E;  %Stores  energy  levels  into  an  array 

53 


a  =  a+1; 

%plottmg  the  wavefunctions - 

hold  on 

plot((X*lel0),((psil/20)+E)) 
plot((Y*lel0),((psi2/20)+E)) 
plot((Z*lel0),((psi3/20)+E)) 
plot((T*  1  e  1 0),((psi4/20)+E)) 

end 

end 

%This  program  finds  the  energy  levels  in  the  Valence  Band 
valenceInGaAsP 

%Drawing  the  finite  well - 

well2  %This  program  draws  the  Conduction  band  well 

% - End  of  aqw  InGaAsP - 
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%This  program  is  titled  Psi  aqw 
%Created:  January  2002 
%Author:  Kevin  Lantz 


%This  is  a  program  that  is  used  to  find  the  wavefunctions  in  each  region  of  the  quantum  well. 
%The  bulk  of  this  program  is  trying  to  find  the  normalization  constant  hi,  which  determines  the 
%normalized  wavefunctions  and  the  coefficients  A1,B1,A2,B2,A3,B3,A4,B4... 

deltax  =  le-10; 

%A11  of  the  coefficients  can  be  written  in  terms  of  B1 

B1  =  1;  %dummy  variable. .. actual  B1  is  Bl*bl  ==  hi 
A2  =  ml2*Bl; 

B2  =  m22*Bl; 

A3  =  (m_ll*ml2  +  m_12*m22)*Bl; 

B3  =  (m_21*ml2  +  m_22*m22)*Bl; 

A4  =  (m_ll_*(m_ll*ml2  +  m_12*m22)  +  m_12_*(m_21*ml2  +  m_22*m22))*Bl; 

%Defming  the  regions  of  the  quantum  well - 

X  =  -(50*deltax):deltax:0; 

Y  =  0:deltax:L2; 

Z  =  L2:deltax:L3; 

T  =  L3:deltax:L3+(50*deltax); 

%Defining  the  unnormalized  wavefunctions - 

psil  =  Bl*exp(-(i*kl*X)); 
psi2  =  A2*exp(i*k2*Y)  +  B2*exp(-(i*k2*Y)); 
psi3  =  A3*exp(i*k3*Z)  +  B3*exp(-(i*k3*Z)); 
psi4  =  A4*exp(i*k4*T); 

%Integrating - 

psil  =  conj(psil); 
psi_2  =  conj(psi2); 
psi_3  =  conj(psi3); 
psi_4  =  conj(psi4); 

%  the  integral  over  all  space  of  |psi|^2  ==  psi*conj(psi)  has  to  be  1 

integrandl  =  (psil.*psi_l)*deltax; 
integrand2  =  (psi2.*psi_2)*deltax; 
integrand3  =  (psi3.*psi_3)*deltax; 
integrand4  =  (psi4.*psi_4)*deltax; 

integral  1  =  sum(integrandl); 
integral2  =  sum(integrand2); 
integraB  =  sum(integrand3); 
integral4  =  sum(integrand4); 

fmali  =  integral  1  +  integraB  +  integraB  +  integraM; 

hi  =  l/(sqrt(fmali));  %normahzation  constant 

%  normalized  wavefunctions  are  just  the  unnormalized  wavefunctions  multiplied  by  the 
normalization  constant 
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Psil  =  bl*psil; 

Psi2  =  bl*psi2; 

Psi3  =  bl*psi3; 

Psi4  =  bl*psi4; 

%  This  part  of  the  program  integrates  the  normalized  wavefunction  in  order  to  check 
%that  it  equals  1 . 

%This  is  used  to  make  sure  that  the  program  is  working... 

Psil  =  conj(Psil); 

Psi_2  =  conj(Psi2); 

Psi_3  =  conj(Psi3); 

Psi_4  =  conj(Psi4); 

Integrandl  =  (Psil.*Psi_l)*deltax; 

Integrand2  =  (Psi2.*Psi_2)*deltax; 

Integrand3  =  (Psi3.*Psi_3)*deltax; 

Integrand4  =  (Psi4.*Psi_4)*deltax; 

Integral  1  =  sum(Integrandl); 

Integral2  =  sum(Integrand2); 

IntegraB  =  sum(Integrand3); 

IntegraM  =  sum(Integrand4); 

fmall  =  Integral  1  +  IntegraB  +  IntegraB  +  IntegraM; 

% - End  of  Psi  aqw - 
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%This  program  is  titled  E  lnGaAsP 
%Created:  February  2002 
%Author:  Kevin  Lantz 

%This  program  has  the  same  form  as  aqw  InGaAsP  except  it  uses  Airy  functions 
clc 

F  =  input('Electric  Field  in  [V/m]:  '); 

Xl_  =  input('Percentage  of  A1  (0-1)  in  region  1 :  '); 

X2_  =  input('Percentage  of  A1  (0-1)  in  region  2:  '); 

X3_  =  input('Percentage  of  A1  (0-1)  in  region  3:  '); 

X4_  =  X1_; 

width  1  =  Inf;  %Extends  from  -Infinity  to  0; 

width2  =  input(' Width  of  the  region  0->a  in  angstroms:  ');  %Well  width 

width3  =  input('Width  of  the  region  a->b  in  angstroms:  '); 

width4  =  Inf;  %Extends  from  L  to  Infinity; 

% - 

y  =  0;  %InitiaI  value  of  y  to  correct  the  logic  of  the  later  condition 
a=  1; 

% - % 

El  =  widthl*le-10;  %Change  input  of  angstroms->meters 
L2  =  width2*le-10; 

L3  =  (width2-t-width3)*le-10; 

L4  =  width4*le-10; 

hbar=  1.055e-34;  %[J*s] 
m0  =  9.11e-31;%[kg] 
e=  1.602e-19;%[C] 

%Conditions  used  to  make  the  program  function  correctly  under  positive  or  negative  bias 
ifF<0; 
zl  =L3; 
z2  =  L2; 
z3  =  0; 

XI  =X4_; 

X2  =  X3_; 

X3  =  X2_; 

X4  =  X1_; 
else 
zl  =0; 
z2  =  L2; 
z3  =  L3; 

XI  =X1_; 

X2  =  X2_; 

X3  =  X3_; 

X4  =  X4_; 
end 

% - % 

ifF<0 

egl  =  1.4236*(1-X1)  0.81612*X1  -  0.13*(X1*(1-X1)); 

eg3  =  1.519*(1-X2)  0.417*X2  -  0.477*(X2*(1-X2)); 

eg2  =  1.4236*(1-X3)  0.81612*X3  -  0.13*(X3*(1-X3)); 
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eg4=  1.4236*(1-X1)  +  0.81612*X1  -  0.13*(X1*(1-X1)); 
m_effl  =  l/((Xl/(0.0365*m0))  +  ((l-Xl)/(0.0795*m0))); 
m_efG  =  l/((X2/(0.026*m0))  +  ((l-X2)/(0.067*m0))); 
m_eff2  =  l/((X3/(0.0365*m0))  +  ((l-X3)/(0.0795*m0))); 
m_eff4  =  l/((Xl/(0.0365*m0))  +  ((l-Xl)/(0.0795*m0))); 


else 

egl  =  1.4236*(1-X1)  +  0.81612*X1  -  0.13*(X1*(1-X1)); 
eg2  =  1.519*(1-X2)  +  0.417*X2  -  0.477*(X2*(1-X2)); 
eg3  =  1.4236*(1-X3)  +  0.81612*X3  -  0.13*(X3*(1-X3)); 
eg4  =  1.4236*(1-X1)  +  0.81612*X1  -  0.13*(X1*(1-X1)); 
m_effl  =  l/((Xl/(0.0365*m0))  +  ((l-Xl)/(0.0795*m0))); 
m_eff2  =  l/((X2/(0.026*m0))  +  ((l-X2)/(0.067*m0))); 
m_eff3  =  l/((X3/(0.0365*m0))  +  ((l-X3)/(0.0795*m0))); 
m_eff4  =  l/((Xl/(0.0365*m0))  +  ((l-Xl)/(0.0795*m0))); 
end 


VI  =0.6*(egl  -0.8161193); 

V2  =  0.6*(eg2  -  0.8161 193);  %Clioose  the  origin  as  the  bottom  left  comer  of  region  2. 
V3  =  0.6*(eg3  -0.8161193); 

V4  =  V1; 

% - % 

Cl  =  (((2*m_effl*e)/((F*hbar)^2))^(l/3)); 

C2  =  (((2*m_eff2*ey((F*hbar)^2))^(l/3)); 

C3  =  (((2*m_efG*ey((F*hbar)^2))^(l/3)); 

C4  =  (((2*m_eff4*e)/((F*hbar)^2))^(l/3)); 

sigmal  =  ((m_effl*C2)/(m_eff2*Cl)); 
sigma2  =  ((m_eff2*C3)/(m_efG*C2)); 
sigma3  =  ((m_efG*C4)/(m_eff4*C3)); 


% - % 

ifF<0; 

DD  =  0; 

V2_  =  V1; 
else 

DD  =  -F*L2; 

V2_  =  V1  -F*L3; 
end 

% - % 

forE=DD:0.001:V2_; 

etal  =  (E  -  VI); 
eta2  =  (E  -  V2); 
eta3  =  (E  -  V3); 
eta4  =  (E  -  V4); 


Alphal  =  -Cl*(F*zl  -I-  etal); 
Alpha2  =  -C2*(F*zl  -I-  eta2); 
Alpha2_  =  -C2*(F*z2  -I-  eta2); 
Alpha3  =  -C3*(F*z2  -I-  eta3); 
Alpha3_  =  -C3*(F*z3  -I-  eta3); 
Alpha4  =  -C4*(F*z3  -I-  eta4); 


%■ 


% 


58 


Ewell2 

o/o - 


mil  =  (pi)*(airy(Alpha2)*airy(3,Alphal)  -  sigmal*(airy(2,Alphal)*airy(l,Alpha2))); 
ml2  =  (pi)*(airy(2,Alpha2)*airy(3,Alphal)  -  sigmal*(airy(2,Alphal)*airy(3,Alpha2))); 
m21  =  (pi)*(sigmal*(airy(Alphal)*airy(l,Alpha2))  -  airy(l,Alphal)*airy(Alpha2)); 
m22  =  (pi)*(sigmal*(airy(Alphal)*airy(3,Alpha2))  -  airy(l,Alphal)*airy(2,Alpha2)); 

m  il  =  (pi)*(airy(Alpha3)*airy(3,Alpha2_)  -  sigma2*(airy(2,Alpha2_)*airy(l,Alpha3))); 
m_12  =  (pi)*(airy(2,Alpha3)*airy(3,Alpha2_)  -  sigma2*(airy(2,Alpha2_)*airy(3,Alpha3))); 
m_21  =  (pi)*(sigma2*(airy(Alpha2_)*airy(l,Alpha3))  -  airy(l,Alpha2_)*airy(Alpha3)); 
m_22  =  (pi)*(sigma2*(airy(Alpha2_)*airy(3,Alpha3))  -  airy(l,Alpha2_)*airy(2,Alpha3)); 

m_ll_  =  (pi)*(airy(Alpha4)*airy(3,Alpha3_)  -  sigma3*(airy(2,Alpha3_)*airy(l,Alpha4))); 
m_12_  =  (pi)*(airy(2,Alpha4)*airy(3,Alpha3_)  -  sigma3*(airy(2,Alpha3_)*airy(3,Alpha4))); 
m_21_  =  (pi)*(sigma3*(airy(Alpha3_)*airy(l,Alpha4))  -  airy(l,Alpha3_)*airy(Alpha4)); 
m_22_  =  (pi)*(sigma3*(airy(Alpha3_)*airy(3,Alpha4))  -  airy(l,Alpha3_)*airy(2,Alpha4)); 

Ml  =[mll,ml2;m21,m22]; 

M2  =  [m_ll,m_12;m_21,m_22]; 

M3  =  [m_ll_,m_12_;m_21_,m_22J; 

M2M3  =  M2*M3; 

M  =  M1*M2*M3; 

M22  =  real(M(2,2)); 

Epsiaqw 

z  =  y; 
y  =  M22; 

ifz.*y  <0 
disp(E) 

Eout(a)  =  E; 
a  =  a  +  1; 
hold  on 

plot((X*  1  e  1 0),((Psi  1  ./(4e5))+E)) 
plot((Y*lelO),((Psi2./(4e5))+E)) 
plot((Z*  1  e  1 0),((Psi3  ./(4e5))+E)) 
plot((T*  1  e  1 0),((Psi4./(4e5))+E)) 


end 

end 

E  valence  InGaAsP 


■End  of  E  InGaAsP- 
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%This  program  is  titled  absorption. 

%Created:  February  2002 
%Author:  Kevin  Lantz 

%This  program  uses  the  wavefunctions  to  find  the  bound-to-bound  and  bound-to-continuum 
%absorption  spectra. 

clc 

aqwIuGaAsP 


% - Psi  at  Eout(l) 


kEl_l  =  sqrt((2*m_effl*(Eout(l)-Vl)*q)/(hbar^2)); 
kEl_2  =  sqrt((2*m_eff2*(Eout(l)-V2)*q)/(hbar^2)); 
kEl_3  =  sqrt((2*m_eff3*(Eout(l)-V3)*q)/(hbar^2)); 
kEl_4  =  sqrt((2*m_eff4*(Eout(l)-V4)*q)/(hbar^2)); 

mEl_ll  =  ((1/2)*((1  -I- ((kEl_l*m_eff2)/(kEl_2*m_effl)))*exp(i*(xl*(kEl_l  -kEl_2))))); 
mEl_12  =  ((1/2)*((1  -  ((kEl_l*m_eff2)/(kEl_2*m_effl)))*exp(-(i*(xl*(kEl_l  -I- kEl_2)))))); 
mEl_21  =  ((1/2)*((1  -  ((kEl_l*m_eff2)/(kEl_2*m_effl)))*exp(i*(xl*(kEl_l  -I- kEl_2))))); 
mEl_22  =  ((1/2)*((1  -I-  ((kEl_l*m_eff2)/(kEl_2*m_effl)))*exp(-(i*(xl*(kEl_l  -  kEl_2)))))); 

m_El_ll  =  ((1/2)*((1  -I-  ((kEl_2*m_eff3)/(kEl_3*m_eff2)))*exp(i*(x2*(kEl_2  -  kEl_3))))); 
m_El_12  =  ((1/2)*((1  -  ((kEl_2*m_efG)/(kEl_3*m_eff2)))*exp(-(i*(x2*(kEl_2  -I- kEl_3)))))); 
m_El_21  =  ((1/2)*((1  -  ((kEl_2*m_efG)/(kEl_3*m_eff2)))*exp(i*(x2*(kEl_2  -I- kEl_3))))); 
m_El_22  =  ((1/2)*((1  -I-  ((kEl_2*m_eff3)/(kEl_3*m_eff2)))*exp(-(i*(x2*(kEl_2  -  kEl_3)))))); 

m_El_l  1_  =  ((1/2)*((1  -I-  ((kEl_3*m_eff4)/(kEl_4*m_eff3)))*exp(i*(x3*(kEl_3  -  kEl_4))))); 
m_El_12_  =  ((1/2)*((1  -  ((kEl_3*m_eff4)/(kEl_4*m_eff3)))*exp(-(i*(x3*(kEl_3  -I- kEl_4)))))); 
m_El_21_  =  ((1/2)*((1  -  ((kEl_3*m_eff4)/(kEl_4*m_eff3)))*exp(i*(x3*(kEl_3  -I- kEl_4))))); 
m_El_22_  =  ((1/2)*((1  -I-  ((kEl_3*m_eff4)/(kEl_4*m_eff3)))*exp(-(i*(x3*(kEl_3  -  kEl_4)))))); 

ME1_1  =  [mEl_ll,mEl_12;mEl_21,mEl_22]; 

ME1_2  =  [m_El_ll,m_El_12;m_El_21,m_El_22]; 

ME1_3  =  [m_El_ll_,m_El_12_;m_El_21_,m_El_22J; 

M2M1  =ME1_2*ME1_1; 

M3M2M1  =ME1_3*ME1_2*ME1_1; 

AE1_1  =  0; 

BE1_1  =  1; 

AE1_2  =  ME1_1(1,1)*AE1_1  ME1_1(1,2)*BE1_1; 

BE1_2  =  ME1_1(2,1)*AE1_1  ME1_1(2,2)*BE1_1; 

AE1_3  =  M2M1(1,1)*AE1_1  M2M1(1,2)*BE1_1; 

BE1_3  =M2M1(2,1)*AE1_1  -HVI2M1(2,2)*BE1_1; 

AE1_4  =  M3M2M1(1,1)*AE1_1  M3M2M1(1,2)*BE1_1; 

BE1_4  =  0; 


deltax  =  le-10; 


X  =  -L3:deltax:0; 

Y  =  0:deltax:L2; 

Z  =  L2:deltax:L3; 

T  =  L3:deltax:(2*L3); 

psiEl  l  =  AEl_l*exp(i*kEl_l*X)  -I- BEl_l*exp(-(i*kEl_l*X)); 
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psiEl_2  =  AEl_2*exp(i*kEl_2*Y)  +  BEl_2*exp(-(i*kEl_2*Y)); 
psiEl_3  =  AEl_3*exp(i*kEl_3*Z)  +  BEl_3*exp(-(i*kEl_3*Z)); 
psiEl_4  =  AEl_4*exp(i*kEl_4*T)  +  BEl_4*exp(-(i*kEl_4*T)); 

psiEll  =  conj(psiEl_l); 
psi_El_2  =  conj(psiEl_2); 
psi_El_3  =  conj(psiEl_3); 
psi_El_4  =  conj(psiEl_4); 

deltalEll  =  (psi_El_l.*psiEl_l)*deltax; 
deltaIEl_2  =  (psi_El_2.*psiEl_2)*deltax; 
deltaIEl_3  =  (psi_El_3.*psiEl_3)*deltax; 
deltaIEl_4  =  (psi_El_4.*psiEl_4)*deltax; 

IE  11  =  sum(deltaIEl_l); 

IE1_2  =  sum(deltaIEl_2); 

IE1_3  =  sum(deltaIEl_3); 

IE1_4  =  sum(deltaIEl_4); 

finalEl  =  IE1_1+IE1_2+IE1_3+IE1_4; 

bEl  =  l/(sqrt(finalEl)); 

PsiEl_l  =bEl*psiEl_l; 

PsiEl_2  =  bEl*psiEl_2; 

PsiEl_3=bEl*psiEl_3; 

PsiEl_4  =  bEl*psiEl_4; 

% - derivative  of  Psi  at  Eout(l) - 

psiEEl  l  =  AEl_l*exp(i*kEl_l*X)  -  BEl_l*exp(-(i*kEl_l*X)); 
psiEEl_2  =  AEl_2*exp(i*kEl_2*Y)  -  BEl_2*expHi*kEl_2*Y)); 
psiEEl_3  =  AEl_3*exp(i*kEl_3*Z)  -  BEl_3*exp(-(i*kEl_3*Z)); 
psiEEl_4  =  AEl_4*exp(i*kEl_4*T)  -  BEl_4*exp(-(i*kEl_4*T)); 

PsiEEll  =bEl*psiEEl_l; 

PsiEEl_2  =  bEl*psiEEl_2; 

PsiEEl_3  =bEl*psiEEl_3; 

PsiEEl_4  =  bEl*psiEEl_4; 

% - Psi  at  Eout(2) - 

kE2_l  =  sqrt((2*m_effl*(Eout(2)-Vl)*q)/(hbar^2)); 
kE2_2  =  sqrt((2*m_eff2*(Eout(2)-V2)*q)/(hbari'2)); 
kE2_3  =  sqrt((2*m_efO*(Eout(2)-V3)*q)/(hbari'2)); 
kE2_4  =  sqrt((2*m_eff4*(Eout(2)-V4)*q)/(hbar^2)); 

mE2_ll  =  ((1/2)*((1  +  ((kE2_l*m_eff2)/(kE2_2*m_effl)))*exp(i*(xl*(kE2_l  -  kE2_2))))); 
mE2_12  =  ((1/2)*((1  -  ((kE2_l*m_eff2)/(kE2_2*m_effl)))*exp(-(i*(xl*(kE2_l  +  kE2_2)))))); 
mE2_21  =  ((1/2)*((1  -  ((kE2_l*m_eff2)/(kE2_2*m_effl)))*exp(i*(xl*(kE2_l  +kE2_2))))); 
mE2_22  =  ((1/2)*((1  +  ((kE2_l*m_eff2)/(kE2_2*m_effl)))*exp(-(i*(xl*(kE2_l  -  kE2_2)))))); 

m_E2_ll  =  ((1/2)*((1  +  ((kE2_2*m_eff3)/(kE2_3*m_eff2)))*exp(i*(x2*(kE2_2  -  kE2_3))))); 
m_E2_12  =  ((1/2)*((1  -  ((kE2_2*m_efG)/(kE2_3*m_eff2)))*exp(-(i*(x2*(kE2_2  +  kE2_3)))))); 
m_E2_21  =  ((1/2)*((1  -  ((kE2_2*m_efG)/(kE2_3*m_eff2)))*exp(i*(x2*(kE2_2  +kE2_3))))); 
m_E2_22  =  ((1/2)*((1  +  ((kE2_2*m_eff3)/(kE2_3*m_eff2)))*exp(-(i*(x2*(kE2_2  -  kE2_3)))))); 

61 


m_E2_l  1_  =  ((1/2)*((1  +  ((kE2_3*m_eff4)/(kE2_4*m_efD)))*exp(i*(x3*(kE2_3  -  kE2_4))))); 
m_E2_12_  =  ((1/2)*((1  -  ((kE2_3*m_eff4)/(kE2_4*m_efD)))*exp(-(i*(x3*(kE2_3  +  kE2_4)))))); 
m_E2_21_  =  ((1/2)*((1  -  ((kE2_3*m_eff4)/(kE2_4*m_efG)))*exp(i*(x3*(kE2_3  +kE2_4))))); 
m_E2_22_  =  ((1/2)*((1  +  ((kE2_3*m_eff4)/(kE2_4*m_efG)))*exp(-(i*(x3*(kE2_3  -  kE2_4)))))); 


ME2_1  =  [mE2_ll,mE2_12;mE2_21,mE2_22]; 

ME2_2  =  [m_E2_ll,m_E2_12;m_E2_21,m_E2_22]; 

ME2_3  =  [m_E2_ll_,m_E2_12_;m_E2_21_,m_E2_22J; 

M2_M1  =ME2_2*ME2_1; 

M3M2_M1  =ME2_3*ME2_2*ME2_1; 

AE2_1  =  0; 

BE2_1  =  1; 

AE2_2  =  ME2_1(1,1)*AE2_1  +  ME2_1(1,2)*BE2_1; 

BE2_2  =  ME2_1(2,1)*AE2_1  +  ME2_1(2,2)*BE2_1; 

AE2_3  =  M2_M1(1,1)*AE2_1  +  M2_M1(1,2)*BE2_1; 

BE2_3  =  M2_M1(2,1)*AE2_1  +  M2_M1(2,2)*BE2_1; 

AE2_4  =  M3M2_M1(1,1)*AE2_1  +  M3M2_M1(1,2)*BE2_1; 
BE2_4  =  0; 

psiE2_l  =  AE2_l*exp(i*kE2_l*X)  +  BE2_l*exp(-(i*kE2_l*X)); 
psiE2_2  =  AE2_2*exp(i*kE2_2*Y)  +  BE2_2*expHi*kE2_2*Y)); 
psiE2_3  =  AE2_3*exp(i*kE2_3*Z)  +  BE2_3*exp(-(i*kE2_3*Z)); 
psiE2_4  =  AE2_4*exp(i*kE2_4*T)  +  BE2_4*expHi*kE2_4*T)); 


psi_E2_l  =  conj(psiE2_l); 
psi_E2_2  =  conj(psiE2_2); 
psi_E2_3  =  conj(psiE2_3); 
psi_E2_4  =  conj(psiE2_4); 


deltaIE2_l  =  (psi_E2_l.*psiE2_l)*deltax; 
deltaIE2_2  =  (psi_E2_2.*psiE2_2)*deltax; 
deltaIE2_3  =  (psi_E2_3.*psiE2_3)*deltax; 
deltaIE2_4  =  (psi_E2_4.*psiE2_4)*deltax; 


IE2_1  =  sum(deltaIE2_l); 

IE2_2  =  sum(deltaIE2_2); 

IE2_3  =  sum(deltaIE2_3); 

IE2_4  =  sum(deltaIE2_4); 

finalE2  =  IE2_1+IE2_2+IE2_3+IE2_4; 

bE2  =  l/(sqrt(finalE2)); 

PsiE2_l  =bE2*psiE2_l; 

PsiE2_2  =  bE2*psiE2_2; 
PsiE2_3=bE2*psiE2_3; 

PsiE2_4  =  bE2*psiE2_4; 

PsiE2_l_  =  conj(PsiE2_l); 

PsiE2_2_  =  conj(PsiE2_2); 

PsiE2_3_  =  conj(PsiE2_3); 

PsiE2_4_  =  conj(PsiE2_4); 
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%■ 


■<Psi(Eout(2))*|k|Psi(Eout(l))>- 


deltaE12_l  =  (PsiE2_l_.*kEl_l.*PsiEEl_l)*deltax; 
deltaE12_2  =  (PsiE2_2_.*kEl_2.*PsiEEl_2)*deltax; 
deltaE12_3  =  (PsiE2_3_.*kEl_3.*PsiEEl_3)*deltax; 
deltaE12_4  =  (PsiE2_4_.*kEl_4.*PsiEEl_4)*deltax; 

E12_l  =  sum(deltaE12_l); 

E12_2  =  sum(deltaE12_2); 

E12_3  =  sum(deltaE12_3); 

E12_4  =  sum(deltaE12_4); 

averageE12  =  E12_l+E12_2+E12_3+E12_4; 

AverageE12  =  averageE12*(conj(averageE12)); 

% - Line  Shape - 

N=  le24; 
eO  =  8.85e-12; 
c  =  3e8; 
n  =  3.5; 

Gamma  =  ((10e-3)*q)/2; 
ap=  1; 

for  EE  =  (Eout(l)-0.1):le-4:(Vl-0.001); 

g  =  (Gamma)/((((Eout(2)-Eout(l)-EE)*q)^2)+((Gamma)^2)); 
gmax  =  (Gamma)/((((0)*q)^2)+((Gamma)^2)); 

C  =  (N*(q^2)*(hbar^3))/((m_effl^2)*e0*n*c); 
alphaE12(ap)  =  (C.*g.*AverageE12)/(EE*q); 

ap  =  ap  +  1 ; 
end 


% 

% - Psi  above  the  well - 

r=l; 

for  Energy  =  (Vl+0.0001):0.0001:(Vl+.025) 


kAl  =  sqrt((2*m_effl*(Energy-Vl)*q)/(hbar''2)); 
kA2  =  sqrt((2*m_eff2*(Energy-V2)*q)/(hbar''2)); 
kA3  =  sqrt((2*m_eff3*(Energy-V3)*q)/(hbar''2)); 
kA4  =  sqrt((2*m_eff4*(Energy-V4)*q)/(hbar''2)); 


mA_ll  =  ((1/2)*((1  +  ((kAl*m_eff2)/(kA2*m_effl)))*exp(i*(xl*(kAl  -  kA2))))); 
mA_12  =  ((1/2)*((1  -  ((kAl*m_eff2)/(kA2*m_effl)))*exp(-(i*(xl*(kAl  +kA2)))))); 
mA_21  =  ((1/2)*((1  -  ((kAl*m_eff2)/(kA2*m_effl)))*exp(i*(xl*(kAl  +kA2))))); 
mA_22  =  ((1/2)*((1  +  ((kAl*m_eff2)/(kA2*m_effl)))*exp(-(i*(xl*(kAl  -  kA2)))))); 

m_A_ll  =  ((1/2)*((1  +  ((kA2*m_efG)/(kA3*m_eff2)))*exp(i*(x2*(kA2  -  kA3))))); 
m_A_12  =  ((1/2)*((1  -  ((kA2*m_efG)/(kA3*m_eff2)))*exp(-(i*(x2*(kA2  +kA3)))))); 
m_A_21  =  ((1/2)*((1  -  ((kA2*m_efG)/(kA3*m_eff2)))*exp(i*(x2*(kA2  +kA3))))); 
m_A_22  =  ((1/2)*((1  +  ((kA2*m_efG)/(kA3*m_eff2)))*exp(-(i*(x2*(kA2  -  kA3)))))); 

m_A_ll_  =  ((1/2)*((1  +  ((kA3*m_eff4)/(kA4*m_efB)))*exp(i*(x3*(kA3  -  kA4))))); 
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m_A_12_  =  ((1/2)*((1  -  ((kA3*m_eff4)/(kA4*m_efB)))*exp(-(i*(x3*(kA3  +kA4)))))); 
m_A_21_  =  ((1/2)*((1  -  ((kA3*m_eff4)/(kA4*m_efB)))*exp(i*(x3*(kA3  +kA4))))); 
m_A_22_  =  ((1/2)*((1  +  ((kA3*m_eff4)/(kA4*m_efB)))*exp(-(i*(x3*(kA3  -  kA4)))))); 

MA_1  =  [mA_ll,mA_12;mA_21,mA_22]; 

MA_2  =  [m_A_ll,m_A_12;m_A_21,m_A_22]; 

MA_3  =  [m_A_ll_,m_A_12_;m_A_21_,m_A_22J; 


M2M1_  =  MA_2*MA_1; 

M3M2M1_  =  MA_3*MA_2*MA_1; 

A1  =  1; 

B1  =  ((-M3M2M1_(2,1))/(M3M2M1_(2,2)))*A1; 
A2  =  MA_1(1,1)*A1  +  MA_1(1,2)*B1; 

B2  =  MA_1(2,1)*A1  +  MA_1(2,2)*B1; 

A3  =  M2M1_(1,1)*A1  +  M2M1_(1,2)*B1; 

B3  =  M2M1_(2,1)*A1  +  M2M1_(2,2)*B1; 

A4  =  M3M2M1_(1,1)*A1  +  M3M2M1_(1,2)*B1; 
B4  =  0; 


PsiAl  =  Al*exp(i*kAl*X)  +  Bl*exp(-(i*kAl*X)); 
PsiA2  =  A2*exp(i*kA2*Y)  +  B2*exp(-(i*kA2*Y)); 
PsiA3  =  A3*exp(i*kA3*Z)  +  B3*exp(-(i*kA3*Z)); 
PsiA4  =  A4*exp(i*kA4*T)  +  B4*exp(-(i*kA4*T)); 

PsiAl_  =  conj(PsiAl); 

PsiA2_  =  conj(PsiA2); 

PsiA3_  =  conj(PsiA3); 

PsiA4_  =  conj(PsiA4); 

% - <Psi(above  well)*  |k|Psi(Eout(  1  ))>- 


deltaEAl  =  (PsiAl_.*kEl_l.*PsiEEl_l)*deltax; 
deltaEA_2  =  (PsiA2_.*kEl_2.*PsiEEl_2)*deltax; 
deltaEA_3  =  (PsiA3_.*kEl_3.*PsiEEl_3)*deltax; 
deltaEA_4  =  (PsiA4_.*kEl_4.*PsiEEl_4)*deltax; 

EA_1  =  sum(deltaEA_l); 

EA_2  =  sum(deltaEA_2); 

EA_3  =  sum(deltaEA_3); 

EA_4  =  sum(deltaEA_4); 

averageEA  =  EA_1+EA_2+EA_3+EA_4; 

AverageEA(r)  =  averageEA*  (conj  (averageEA)); 

% - Finding  the  absorption  coefficient- 

C  =  (N*(q^2)*(hbar^2))/(2*(m_effB2)*eO*n*c); 

CC  =  (((((Energy-Vl))*q)^2)  +  (Gamma^2)); 

rho  =  ((l/CC)^(l/4)); 

C_  =  (((Energy-Eout(l))*q)^2); 
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A_  =  (sqrt((2*m_effl)/C_))*rho; 

alpha(r)  =  C  *  A*  Average  EA(r); 

r  =  r+  1; 
end 

% - End  of  absorption- 
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